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ISOTOPIC COMPOSITION OF SILVER IN AN IRON METEORITE* 


V. Rama Murthy 
School of Science and Engineering, University of California, La Jolla, California 
* (Received November 21, 1960) 


The recent discovery of radiogenic Xe**® in 
Richardton chondrite by Reynolds’ has empha- 
sized the importance of a search for extinct ra- 
dioactivities. Several such possible nuclides 
have been suggested.”** As part of a program of 
investigating the isotopic composition and abun- 
dances of several elements in meteorites, the 
isotopic composition of silver in the troilite of 
Toluca iron meteorite was determined and this 
paper discusses the results. 

Pd'” is an extinct radioactive nuclide which 
decays to Ag’. The half-life of this decay is 
reported to be 7.510® years by Sullivan.* Dur- 
ing the formation and separation of silicate and 
metal phases of meteorites, Pd will be incor- 
porated into the metal phase because of its sider - 
ophile nature. Goldschmidt’ classified silver 
as both a siderophile and chalcophile though the 
observational data for the meteorites are unre- 
liable. Thermodynamic calculations assuming 
ideal solution laws for Ag,S in troilite and silver 
in iron-nickel substantiate a siderophile classifi- 
cation. However, ideal solution laws may not 
apply because elemental silver forms a two- 
liquid system with molten iron and hence the cal- 
culations must indicate a more siderophile char- 
acter than is actually the case. Silver, therefore, 
is likely to be fractionated from Pd in the metal 
phase. If the formation of the metal phase of 
meteorites happened within a few half-lives of 
Pd'” after the end of nucleosynthesis, because 
of the above-mentioned fractionation between Pd 
and Ag, the metallic phase is likely to show an 
enrichment of Ag’” as a result of decay of Pd'”. 
During the subsequent cooling of the metal phase 


and formation of the troilite nodules, it is pos- 
sible that the Ag’”-enriched:silver from the me- 
tal phase was concentrated in the troilite. 

Silver was extracted from the troilite of Toluca 
meteorite by ion exchange techniques (to be pub- 
lished elsewhere). A few micrograms of Ag were 
isolated and mounted on a tantalum filament as 
Ag,S and a drop of 2% NH,NO, was added to the 
filament. The mass spectrometer used was a 12- 
in. radius and 60°-deflection instrument built by 
Nuclide Analysis Associates. 

The results of measurement of the isotopic 
composition of silver in Toluca troilite and of 
silver from reagent grade Baker and Adams 
AgNO, are given in Table I. The data are re- 
ported as measured in both cases without any 
corrections. It is possible th.t our mass spec- 
trometer has a systematic mass discrimination 
in the range of silver isotopes, but this will not 
affect a comparison of terrestrial silver with 


Table I. Isotopic compositions of silver in Toluca 
troilite and terrestrial silver. 








No. of 
Sample Ag!®"/Ag!® ratios 
Toluca troilite 
Run I 1.097 + 0.002 14 
Run II 1.100 + 0.002 15 
Run Il 1.094 + 0.002 20 
Average 1.097 +0. 002 
Reagent AgNO, -I 1.076 + 0.002 22 
Reagent AgNO, - II 1.072 + 0.002 18 
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meteoritic silver; any systematic discrimina- 
tion would affect both measurements equally. 

The isotopic composition of Ag in Toluca troi- 
lite was measured as three sets of data, each set 
comprising 14 to 20 ratios of Ag’”/Ag’™®, during 
a six-hour run of stable emission. The mean 
value of Ag?*/Ag*®® in Toluca troilite is 1.097 
+0.002. Our measurement of Ag’”/Ag’™® ratio of 
terrestrial Ag, using reagent grade AgNO,, is 
within experimental error the same as the value 
of 1.0756 for the Ag standard of the National Bu- 
reau of Standards and several other samples of 
terrestrial Ag as carefully measured by Shields, 
Craig, and Dibeler.® 

Hess, Marshall, and Urey’ analyzed the com- 
position of silver in Toluca troilite and found that 
silver in this meteorite was very similar to 
terrestrial silver. However, their estimated ex- 
perimental error was large enough to mask the 
2% enrichment of Ag’”/Ag’™ ratio in Toluca 
troilite, reported in this work. 

The excess Ag’” in Toluca troilite may have 
been caused by selective fractionation of the light 
isotope of Ag into the troilite during long period 
of cooling of the iron phase, estimated by Urey® 
to be of the order of 10” - 10° years, or, by some 
nuclear processes during the last stages of nu- 
cleosynthesis, in the solar system.® However, 
if the excess Ag’” in the Toluca troilite is attri- 
buted to the decay of Pd’”, the time interval, 

At, between the end of nucleosynthesis and the 
formation of iron meteorites can be estimated. 

In the “sudden synthesis” model, which postu- 
lates element formation in one event, A? is given 
by 

At= Tio7[In(Pd??/Ag**"*) + In(Pd?™ / Pd?) | ‘ 


where 7,,, is the mean life of Pd’, (Pd*°*/Ag**’*) 
is the ratio of Pd'® to radiogenic Ag’™ in the 
meteorite, and (Pd’”/Pd'®) is the ratio of these 
isotopes at the end of nucleosynthesis. The Pd 
content of Toluca iron meteorite is reported to 
be 4.72 ppm by Goldberg, Uchiyama, and Brown. 
Judging from the cosmic abundance curve of 
Suess and Urey" at Pd mass region, the ratio of 
the abundances of two neighboring odd- mass nu- 
clides, Pd'*’/Pd’®, can be taken to be about unity. 
The silver content of iron meteorites is not well 
known. Hess, Marshall, and Urey report less 
than 0.2 ppm in the Canyon Diablo iron and 
Lovering”’ reports less than 0.5 ppm in three 
troilites including Toluca. Assuming a silver 
content of 0.1 ppm in Toluca iron, Af in the 
“sudden synthesis” model is about 8.410" years. 
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In the more commonly accepted “continuous 
synthesis” model of Burbidge, Burbidge, Fowler, 
and Hoyle,** the Pd*” abundance at the end of 
nucleosynthesis is expected to be in steady-state 
proportions. In this model, Aft is given by 


At = 7,9, [In(Pd?*/Ag’*) + In(Pd! / Pd’) 
+ 1n(7,9,/T)}. 


The value of A¢ in this model is 10’ years. This 
estimate of A? is perhaps not very significant be- 
cause the half-life of Pd'” is short compared to 
the total time span of nucleosynthesis, T, of the 
order of 10'° years. The abundance of Pd’” at 
the end of nucleosynthesis will depend to a large 
extent only on the production of this isotope in 
the last supernova. On this basis, assuming a 
Pd’ /Pd'® ratio of 0.1 at the end of a “contin- 
uous model” nucleosynthesis, Af is about 4.5 
x10’ years. This estimate of At is similar to 
the value of 4x10" years obtained for Canyon 
Diablo iron by Anders and Stevens using the 
Pb?-T1?°> decay scheme. These estimates of 
At are still not precise enough to evaluate the 
role of many other shorter lived radioactivities 
in the early development of meteorites and other 
planets. 

Further work is in progress to look for this ef- 
fect in other iron meteorites, particularly those 
which are rich in Pd. 

I wish to thank Professor H. C. Urey, Univer- 
sity of California, La Jolla, for helpful discus- 
sions. 
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HELIUM II IN ROTATION* 


John D. Reppy, David Depatie, and C. T. Lane 
Department of Physics, Yale University, New Haven, Connecticut 
(Received November 21, 1960) 


In a recent experiment,’ Pellam has reported 
some measurements concerned with the velocity 
field of He II when contained in a simply connected 
vessel subject to steady rotation about its verti- 
cal axis of symmetry. These results are inter- 
esting, and unexpected, since they are contrary 
to the predictions of any current theory of liquid 
helium. We felt it important, therefore, to ex- 
amine the problem using a different approach. 

In our experiment the helium is contained ina 
light cylindrical glass vessel (bucket) which can 
be given an angular velocity about the vertical 
axis. The volume occupied by the helium was ap- 
proximately 2.5 cm in diameter by 11 cm in depth. 
The bucket and its contents are suspended from a 
Beams type magnetic bearing, the whole system 
being enclosed in a high vacuum. Under these 
circumstances the parasitic rotational drag on 
the system was less than 5x10~° dyne cm which, 
for our purposes, was negligible. 

The procedure consists in giving the system, 
initially at rest, a sudden (~1/2 sec) rotational 
impulse and thence allowing it to coast freely. 
Observations are then made of the angular velo- 
city of the system as a function of time and Fig. 1 
shows two typical results, one above and the other 
below the A point. As seen, the velocity of the 
system decreases with time and eventually reaches 
a steady value. Since angular momentum is con- 
served, the slow-down of the system must mean 
that the liquid, initially at rest, is slowly acceler- 
ated into rotation by the walls of the bucket. The 
“reaction time” is sharply temperature dependent, 
being the faster the higher the temperature, and 
this indicates that the mechanism of interaction 
between bucket and liquid is different for the nor- 
mal and superfluid components of the helium. 

Control experiments show that no appreciable 
angular momentum is transmitted to the liquid in 
the very short initial period during which the sys- 
tem is being brought into rotation. Hence a sim- 
ple calculation, making use of readily measured 
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FIG. 1. The angular velocity of the system, subject 
to no external torques, as a function of time. Two runs 
are shown, for He I and He Il, respectively, illustra- 
ting the very different types of interaction in the two 
fluids. 


constants, permits us to find the angular momen- 
tum acquired by the helium in reaching its equi- 
librium state. 

The result, in all cases, shows that this angu- 
lar momentum is always classical to within an 
average experimental error of about five percent. 
By classical is meant that, at equilibrium, the 
helium is rotating in solid body type motion with 
the same angular velocity which the bucket then 
possesses~—a situation which would also occur in 
the case of an ordinary viscous fluid. This be- 
havior, of course, is precisely that predicted for 
a pure superfluid by the Onsager-Feynman’ vor - 
tex line theory and this theory has met with some 
success in accounting for a variety of other super- 
fluid experiments. 

Figure 2 is a plot of the ratio of the angular mo- 
mentum so measured (L) to the classical value 
(L;) as a function of temperature. The above 
ratio is identical to Pellam’s pe¢¢/p and his re- 
sults are therefore included on the same graph 
for comparison. 
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compared to the classical value (L,) as a function of 
temperature. Circles are the present experiment and 
the dashed curve is Pellam’s result (reference 1}. The 
rotational velocity in all runs was near 2 rpm. 


The most obvious difference between our tech- 
nique and that of Pellam lies in the fact that his 
Rayleigh disk provides a semimicroscopic meas- 
ure of the velocity field whereas ours provides an 
average over such a field. But this fact alone can 
hardly explain the vast difference in results which 
the two experiments show. It must be pointed out 
that earlier experiments on angular momentum by 


Hall* and Walmsley and Lane* are at variance with 
both of the above experiments and, of course, with 
the theory. Plausible, though not entirely certain, 
reasons can be advanced to explain this discrep- 
ancy with our present result but not, at all obvi- 
ously, with Pellam’s. 

The present experiment is the only one of those 
cited which is in concordance with current theories 
of He II. Further, our method appears to subject 
the helium to much less extraneous disturbance, 
such as induced turbulence in the liquid, than do 
any of the others. We suggest, in consequence, 
that Pellam’s result is possibly connected with 
some such extraneous effect. 





*Assisted by a grant from the National Science Founda- 
tion, 
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Physics, edited by J. C. Gorter (Interscience Publishers, 
New York, 1955), Vol. 1. 

3H. E. Hall, Trans. Roy. Soc. (London) 250, 359 
(1957). 

‘R, H. Walmsley and C, T. Lane, Phys. Rev. 112, 
1041 (1958). 








MAGNETIC MOMENT OF TRANSITION METAL ATOMS IN DILUTE SOLUTION AND THEIR 
EFFECT ON SUPERCONDUCTING TRANSITION TEMPERATURE 


B. T. Matthias, M. Peter, H. J. Williams, A. M. Clogston, E. Corenzwit, and R. C. Sherwood 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received November 15, 1960) 


We have found a decisive experiment which 
shows the conditions under which the magnetic 
moment of a transition metal atom remains lo- 
calized when in dilute solution in another transi- 
tion metal. We have correlated this effect with 
a drastic lowering of the superconducting transi- 
tion temperature in a particular superconductor. 

Some time ago a series of observations was 
reported which showed that the superconducting 
transition temperature of various metals and 
alloys could be greatly decreased by the presence 
of small amounts of rare earth metals.' This ef- 
fect has been identified as due to an exchange 
interaction between the conduction electrons and 
the localized f-shell electrons.” In contrast to 
this behavior, it has been reported that the addi- 
tion of transition metal impurities raises the 
superconducting temperature of titanium by an 
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order of magnitude® and that of zirconium some- 
what less.‘ In the case of the fifth column ele- 
ments vanadium and niobium, dilute transition- 
metal impurities lower the transition tempera- 
ture in a way strictly in accord with the electron 
concentration.® This behavior is consistent with 
the fact that even large concentrations of iron 
dissolved in Ti, V, or Nb do not show any mag- 
netic interaction or temperature-dependent sus- 
ceptibility. As far as we are aware no localized 
moment has been observed for a transition-metal 
impurity in dilute solution in another transition 
metal of the fourth or fifth column. 

In this paper we now report a strong lowering 
of the transition temperature of the superconduct- 
ing alloy Mo,,,Re,,, by the solution of small 
amounts of metals belonging to the first- row 
transition series. These experiments are re- 
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conciled to those previously reported by showing 
that the 3d metals in dilute solution in metals of 
the fourth row acquire a magnetic moment ab- 
ruptly between the fifth and sixth columns. This 
result is consistent with the fact that very dilute 
solutions of iron in palladium are known to be- 
come ferromagnetic with an effective moment 
for the iron exceeding two Bohr magnetons per 
atom.® 

The sixth column element Mo is not a super- 
conductor. To study the effects of 3d impurities 
upon transition metals in the sixth column and 
beyond, we have therefore made use of the 
molybdenum-rhenium alloys whose high super- 
conducting transition temperatures were dis- 
covered by Hulm.’ This temperature decreases 
linearly with iron concentration at a rate of 22° 
per atomic percent Fe. In these concentrations 
the addition of nonmagnetic elements has little 
effect. The specific depression is about 150 
times greater in Mo, ,Re, , than in the neighbor- 
ing elements Nb and V. 

Susceptibility measurements have been made 
on the Mo, ,Re, , alloys with iron as an impurity 
which show a Curie-Weiss law for the suscepti- 
bility.° The constant c in the relation x =c/(T-T,) 
is plotted in Fig. 1 as a function of concentration. 
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FIG. 1. The constant c in y=c/(T -T,) corresponds 
to 2.3ug per atom of iron in MogsRep,,. Tec is nega- 
tive and a smooth function of concentration (broken 
line), 


x x 10° 


From it we derive p = g|S(S + 1)]*=3.2, which for 
S=1 leads to a saturation moment of 2.3 B per 
iron. We conclude that the lowering of the tran- 
sition temperature of Mo, ,Re,,, by these impuri- 
ties is due to an exchange interaction with the 
localized spin and is analogous to the effects - 
produced by rare earth ions in lanthanum and 
other metals and alloys.’ We have also plotted 
in Fig. 1 the Curie temperature 7, which turns 
out to be negative and a smooth function of the 
concentration. 

Both niobium and molybdenum crystallize ina 
body-centered cubic structure and form a con- 
tinuous series of solid solutions. The presence 
of a localized moment for iron atoms dissolved 
in Mo, ,Re, ., and the absence of such a moment 
in Nb, suggests that a localized moment will ap- 
pear somewhere in the series between Nb and 
Mo. In Fig. 2 we show the effective magnetic 
moment per iron atom, ~, as a function of Mo 
concentration, obtained from susceptibility data 
taken on Nb - Mo solid solutions containing 1 
atomic percent iron. An appreciable moment 
appears first at 50% Mo. The Nb - Mo alloys up 
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FIG. 2. Effective magnetic moment per iron atom, 
p, in Nb-Mo solutions containing 1 atomic percent Fe. 
The broken line gives the susceptibility at 300° in the 
absence of Fe. From it is derived the energy scale 
explained in the text. 
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to 30% Nb remain superconducting and no appre- 
ciable reduction in T,; due to the iron was ob- 
served in the present experiments, consistent 
with the conclusion advanced above. 

We have also plotted in Fig. 2 the susceptibility 
X per gram of the Nb - Mo alloys without iron. 

If we assume that the susceptibility measures 
the density of states at the Fermi surface, we 
may relate an increase in the average number of 
electrons per atom N to an increase in the Fermi 
energy Ep by the relation AEr =. p*AN/x M,A, 
where M, is the mass of an atom of unit atomic 
weight and A is the average atomic weight. The 
energy scale derived in this way is shown in 

Fig. 2. 

The sharp rise in magnetic moment per atom 
seen in Fig. 2 is a remarkable new phenomenon 
and provides some insight into the electronic 
structure of iron present as an impurity in the 
Nb - Mo crystal. A possible interpretation of the 
observations is the following. Each iron atom 
has 5d states. Following suggestions by Mott 
and Stevens® that apply to the bec structure, we 
may assume that 3 of these states of symmetry 
to, become part of the conduction band of the 
host lattice. The remaining two states of sym- 
metry e., we suppose to be only weakly coupled 
to surrounding wave functions and therefore es- 
sentially localized. For the alloys of less than 
50% Mo we assume that the Fermi level lies be- 
low the eg states and no moment is observed. 
Above 80% Mo we suppose that the Fermi level 


has crossed these states and that they are occu- 
pied by two electrons whose spins align accord- 
ing to Hund’s rule to give the observed moment. 
The antiparallel spin states will be split off by 
exchange and correlation energy so that they lie 
above the Fermi surface and are unoccupied. 
From the energy scale in Fig. 2 we may esti- 
mate the effective width of the e g band to be 
about 0.3 ev. 
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EFFECT OF SPIN-ORBIT SPLITTING ON THE FERMI SURFACES OF THE 
HEXAGONAL-CLOSE- PACKED METALS* 


Morrel H. Cohen and L. M. Falicov 
Institute for the Study of Metals, University of Chicago, Chicago, Illinois 
(Received November 14, 1960) 


All theories of the electronic properties of 
metals and alloys having the hexagonal-close- 
packed (hcp) structure have supposed that the en- 
ergy gaps vanish across the hexagonal face of 
the first Brillouin zone.’»? The simple argument 
based on the vanishing of the structure factor put 
forward first by Jones*® was not correct for a 
general potential. A rigorous basis for the stick- 
ing together of the bands was provided by the 
symmetry properties and time-reversal argu- 
ments of Herring.‘ These were generalized to 
include spin-orbit coupling by Elliott for various 
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crystal structures.® In particular, one sees 
from his results that the degeneracy is removed 
at the point H and along the symmetry lines S and 
S’, but not at A and L nor along the line R 
(Fig. 1). The implications of Elliott’s work for 
the electronic properties of hcp metals and alloys 
appear to have been overlooked. Two questions 
arise: How big is the splitting, and what modifi- 
cation of the existing theories is required? 

We have extended Elliott’s analysis to a general 
point in the plane, and have found, as expected, 
the degeneracy to be lifted. Near A the splitting 
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FIG. 1. The first Brillouin zone for hcp structures, 
showing points and lines of symmetry in the hexagonal 
face. 


is very small because it is required to vanish by 
symmetry along the three lines R intersecting at 
A. The splitting increases from A along S ini- 
tially as (Ak)* towards a maximum value for the 
entire plane at H. 

We have calculated numerical values of the 
splittings at H for magnesium using the results 
of a previous orthogonal plane wave (OPW) cal- 
culation of the band structure.® A six-OPW ap- 
proximation to the wave function was used; 
fewer OPW’s would not have yielded the correct 
order of magnitude, more would have required 
an inappropriate amount of computation. The 
results are listed in Table I. 

From the detailed wave functions one would 
have anticipated the splitting to be small for H,, 
and for H, and H, to be somewhat greater than 
half the spin-orbit-coupling parameter for a 3p 
electron in atomic Mg (5x10~° ev from the term 
values’ of the configuration 3s3p *P.7). Both are 
found to be the case. Approximately the same 


Table I. The energies at H in magnesium. 








Without spin With spin 





Level Energy (ev) Levels Splitting (ev) 
H, 9.33 ie 2.83 x 107° 
Hs 
H, - 
H, 9.63 s 2.86 x 107° 
Hy 
Hs; 12.22 = ies, 3.1x 10> 


Fermi level at about Ey, =9.25 ev. 














(b) 


FIG. 2. The change in the connectivity properties 
of the Fermi surface of magnesium due to spin-orbit 
coupling. Diagram (a) shows the “without spin” case 
and (b) the surface in the second band when spin-orbit 
coupling is taken into account. 


ratios of splittings in the crystal to atomic spin- 
orbit-coupling parameter should hold for the 
other hcp metals: Be, Zn, Cd, and Tl. For 
these metals, therefore, the spin-orbit splittings 
are 0.05, 9.5, 28, and 129 times as large as the 
corresponding splittings for Mg. 

No major changes of existing alloy theories 
are required because the spin-orbit splitting 
vanishes at the center of the hexagonal face and 
is very small in a sizable region around the cen- 
ter. However, the topological features of the 
Fermi surfaces of the pure metals are funda- 
mentally changed, and the double zone scheme 
ordinarily used in their representation is no 
longer valid. Without spin-orbit coupling, that 
piece of the Fermi surface corresponding to 
holes in the first and second bands in Mg® and 
Zn’ is multiply connected with infinite extent 
normal to the c axis but not parallel to it 
[Fig. 2(a)]. With such a Fermi surface the ob- 
served lack of saturation of the magnetoresist- 
ance for fields normal to the c axis®»® could not 
be understood. The spin-orbit splitting, however, 
separates the small pockets of holes in the first 
band from the region of holes in the second and 
causes the latter to extend infinitely along the c 
axis [Fig. 2(b)]. This explains the observed 
magnetoresistance, because now open orbits 
exist for all directions of magnetic field normal 
to the hexad axis. 
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VOIGT EFFECT IN SEMICONDUCTORS 


S. Teitler and E. D. Palik 
U. S. Naval Research Laboratory, Washington, D. C. 
(Received October 10, 1960) 


Double refraction of an electromagnetic wave 
passing through an isotropic substance under the 
influence of a static magnetic field in a direction 
normal to the direction of propagation is known 
as the Voigt effect." We take as the equation of 
motion for an electron in the conduction band, 

at 


dt + dr = 
ga 10/1) $= se/m*(E + xii), (1) 





2n * = € - 40,7 /B[(a/8)* + w?7?] 


where fF is the displacement vector of the electron, 
T is the mean collision time of the electron, E is 
the incident radiation field, H is the external mag- 
netic field, and m* is the effective mass of the 
electron. 

A straightforward classical calculation using 
Maxwell’s equations yields, using unrationalized 
cgs units, 





_ 4n0,T ’ 1610,"a? = 2 
{( [azar so) * OB (a /B) + o"7"P 
2K ° = -€ +470,7/B[(a/B)? + wt? ] 

4no,T 16770,?a? (2b) 


1 v2 
+{( “aera pe 


where the magnetic permeability has been assumed 
equal to one, n, is the index of refraction for radi- 
ation with electric vector normal to both the mag- 
netic field and the direction of propagation, K, is 
the extinction coefficient in the normal direction, 
and 


a a1 {uty 2r2/lut+(u* -0, 20" 
‘ 2 
B=1 “-. (w? - wr Ae +(w* - wr" lh, 


a, =Ne?1/m*, o.* eH /m*c, w," =4nNe?/m*e, 


0 


with N equal to the free electron concentration. 
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For radiation with electric vector parallel to the 
magnetic field, the index of refraction, nj, and 
the extinction coefficient, K\, follow immediately 
from Eqs. (2) by setting a =8 =1. 

The phase shift between components parallel 
and normal to the magnetic field when incident 
radiation of wavelength A traverses a substance 
of thickness d in a direction normal to the mag- 
netic field is 


6=2nd(n \ -n )/r. (3) 


This phase shift takes a particularly simple form 
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when 
wr>>1, w>w, ww, 
c p 
is. 
i Nde*H?,$ w? || w? | (4) 
es ste So oP 
4n*e”*m'"*c” |w 20, al ¥, w. 


where we have taken (w,”/w”) and w,”/w* to be 
small quantities of the same order of magnitude, 
X. Equation (4) is correct to terms of order X°. 
To check Eq. (4) experimentally we used a sam- 
ple of n-type InAs with carrier concentration of 
10'7/cm’ which is thereby extrinsic at room tem- 
perature. The mobility was in excess of 2 x10* 
cm?/volt-sec. Under experimental conditions 
employed, the validity of Eq. (4) is assured. Also 
since wt >>1, this sample had the added conven- 
ience that away from the cyclotron resonance 
region, there was essentially no difference be- 
tween K, and K\ to complicate the double-refrac- 
tion experiment. Infrared radiation at various 
wavelengths in the region 14u- 19y and plane 
polarized at an angle of 45° to the magnetic field 
with a pile-of-plates polyethylene polarizer was 
incident on the InAs at room temperature ina 
direction normal to the magnetic field. We de- 
tected circularly polarized radiation with a CsI 
Fresnel rhomb and a second polyethylene polari- 
zer as a function of magnetic field. Fields as 
high as 70 kilogauss were used. A typical result 
is shown in Fig. 1. The first peak, which signifies 
a phase shift of 7/2, occurred in all cases before 
the cyclotron resonance absorption began and may 
be used directly in Eq. (4). The same results 
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FIG. 1. Transmission of InAs versus magnetic field 
at 17.9 microns, illustrating Voigt effect. 


were obtained when the direction of the varying 
field was reversed, thereby indicating that the 
double refraction is an even function of the mag- 
netic field. While the occurrence of extrema in 
the observation of circularly polarized radiation is 
unaffected by any difference in extinction coeffi- 
cients, a correction for the changing background 
due to cyclotron resonance absorption must be 
made. In the raw data as shown in Fig. 1, this 
cyclotron resonance absorption makes it difficult 
to assign a reliable magnetic field to the minimum. 

In order to establish the properties of the double 
refraction we restricted ourselves to the use of 
the first maxima which are essentially unaffected 
by cyclotron resonance absorption. The effective- 
mass ratio obtained by inserting these data into 
Eq. (4) is plotted in Fig. 2, as a function of mag- 
netic field. The intercept at 0.0313 for zero mag- 
netic field is consistent with the low-field Faraday- 
rotation effective-mass ratio obtained for the 
same sample. This compares with an effective- 
mass ratio of 0.023 obtained on pure InAs by cy - 
clotron resonance.” This indicates that we are 
measuring an effective mass substantially above 
the bottom of the nonparabolic energy band, which 
is to be expected since our sample is relatively 
impure. This is consistent with Austin’s*® result 
of 0.043 for the effective-mass ratio in a more 
impure sample of InAs with a carrier concentra- 
tion of 1.04x10"*/cm*. The slope of the effective- 
mass ratio curve of Fig. 2 is also consistent with 
the slope of the effective-mass ratio versus mag- 
netic field for cyclotron resonance data for pure 
material.” 

In preliminary experiments on InSb we have 
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FIG. 2. Variation of effective-mass ratio with mag- 
netic field for InAs at room temperature. 


547 








VOLUME 5, NUMBER 12 


PHYSICAL REVIEW LETTERS 





DECEMBER 15, 1960 





also obtained the Voigt effect using the above 
procedure. However, for more general appli- 
cation to semiconductors with heavier masses, 
it is desirable to modify the procedure. This is 
so because the phase shifts are less than 7/2 for 
wavelengths for which the materials are trans- 
parent and for the available static magnetic fields. 
We again start with radiation polarized 45° to the 
magnetic field but now observe the major axis, 
the ratio of axes, and the orientation of the re- 
sultant elliptically polarized radiation due to 
double refraction for a fixed magnetic field and 
wavelength by rotating a linear polarizer before 
the detector. 

Recently Gurevic et al.* have independently dis- 
cussed magnetic double refraction properties of 


semiconductors which may be similar to those 
considered in the present Letter. 

We gratefully acknowledge many informative 
discussions with Dr. R. F. Wallis of this Labo- 
ratory. We also acknowledge an interesting dis- 
cussion concerning magneto-optic effects with 
Dr. R. A. Toupin. j 





1See, e. g., H. A. Lorentz, The Theory of Electrons 
(G. E. Stechert Company, New York, 1916), 2nd ed., 
p. 164. 

2E. D. Palik and R. F. Wallis (to be published). 

31. G. Austin, J. Electronics and Control 8, 167 
(1960). ) 

‘L. D. Gurevié, I. P. Ipatova, and Z. I. Uricky, 
International Conference on Semiconductor Physics, 
Prague, 1960 (unpublished). 








EFFECT OF LANDAU LEVELS UPON TUNNEL CURRENTS IN INDIUM ANTIMONIDE 


A. G. Chynoweth, R. A. Logan, and P. A. Wolff 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received November 7, 1960) 


It has already been shown’~* that when a mag- 
netic field, H, is applied to an InSb Esaki junc- 
tion there is a reduction in the tunnel current 
due to the change in the tunnelling probability 
that arises primarily through the increase in the 
energy gap with H. An additional effect has now 
been observed at low temperatures, namely, an 
oscillation of the tunnel current with H. 

The junctions used in this work were made by 
alloying pellets of 0.1% Cd in In into n-type InSb 
containing an average excess donor density of 
10’ cm~*. The excess acceptor concentration on 
the p side was estimated to be 10'* cm™*. Re- 
sults obtained from two junctions are described 
in this paper; one was an Esaki junction, the 
other a backwards diode. Although these junc- 
tions were made from the same crystal, it is 
believed that the well-known inhomogeneities in 
the donor concentration were responsible for the 
observed spread in the junction characteristics. 
Results basically similar to those described here 
were obtained with several other junctions. 

Current-voltage characteristics taken at 4°K 
showed that as H was increased from zero to 
85000 gauss, the tunnel current underwent an 
oscillation superimposed on the decrease of cur- 
rent noted previously. The oscillations were ob- 
served with H||E, the junction field, and not with 
H1E. For more detailed studies, a constant 
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current was passed through the junction and the 
potential drop across the junction was plotted 
continuously as a function of H on an X-Y re- 
corder. (Plotting voltage was virtually equiva- 
lent to plotting current since the current was 
roughly proportional to the bias for small 
changes in the latter.) Figure 1 shows the re- 
sults obtained with the Esaki junction at small 
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FIG. 1. The bias required to maintain a constant 
tunnel current as a function of magnetic field for two 
different junctions: Junction E was an Esaki diode; 
junction B was a backwards diode. The oscillations in 
bias are attributed to the Landau levels in the conduc- 
tion band. 
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forward bias and with the backwards diode at 
larger biases, the forward bias being such as to 
operate the junction in the excess current range. 
It is apparent that as H increases, small vol- 
tage oscillations appear, that the amplitudes and 
field spacings of these oscillations increase with 
H, and that the oscillations are similar for both 
forward and reverse biases. 

To establish whether the oscillations were 
caused by a series magnetoresistance effect, the 
location of the phonon-bump on the current-vol- 
tage characteristics was examined as a function 
of H using equipment* that automatically plotted 
(d*I/dV?) against V. At H=0 the phonon bump 
occurred at a forward bias of 24.54 0.5 mv, in 
agreement with other measurements,’ and it re- 
mained at this value to within 1.0 mv up to the 
highest values of H, thereby eliminating series 
resistance effects as the cause of the voltage 
fluctuations. This is in line with previous ob- 
servations.’ 

It is possible to ascribe the oscillations to the 
oscillatory variation in the density of electrons 
available for tunnelling brought about by the 
Landau quantization in the conduction band on the 
n-type side. At low biases most of the tunnelling 
transitions are to and from energy levels located 
near the Fermi levels on both sides of the junc- 
tion. For electrons and light holes the density 
of states at the Fermi surface is proportional to 


(2e,+€0) 


2, fe -+ EO)Er - Eglin +2)hw + ugH|}?’ 


(1) 





where €y is the Fermi energy, €¢ the energy 
gap, w, the cyclotron frequency, and yg the spin 
moment. This function has a sharp maximum 
whenever the bottom of a Landau sub-band 
crosses the Fermi level and thereby gives rise 
to fluctuations in the tunnel current. 

Expression (1) predicts that if H~* for the nth 
maximum is plotted against , a straight line 
should result. Such plots are shown for the two 
junctions in Fig. 2 where it is apparent that very 
good straight lines are obtained. The slope ofthe 
line, which is proportional to m*ep(er + €g)/Eg, 
is greater for the Esaki junction than for the 
backwards junction by a factor of 2.4. This ratio 
can also be estimated from the current densities 
of the two diodes. At small biases, the current 
density in the Esaki junction was 70 times that 
of the backwards diode so that the junction width 
constants differed by about 470 A whereas their 
average width was about 1440 A. Hence, the 
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FIG. 2. Plots of the oscillation number versus H™ 
derived from the curves of Fig. 1. 


donor concentration for the Esaki junction was 
about 2.5 times that for the backward junction, 
thereby making the predicted ratio of the slopes 
about 1.9, in reasonable agreement with the ob- 
served value. The values of m* derived from the 
two lines of Fig. 2 were equal to within +5%, 
their average value of 0.011m, being in reason- 
able agreement with other determinations of the 
electron mass at the bottom of the conduction 
band,*® considering the uncertainties in the esti- 
mates of €pr. 

These data show that the oscillations in the 
tunnel current are due to fluctuations in level 
density of the Fermi surface in the n-type sides 
of the diodes only. It is important to realize that 
these current maxima [see Eq. (1)] are due to 
carriers whose energy is entirely transverse to 
the junction field. Thus the observation of oscilla- 
tions indicates that an appreciable fraction of the 
tunnel current is carried by electrons whose 
velocity is, initially, perpendicular to the field. 
In electron-light hole tunnelling, this transverse 
energy must be conserved. From the known band 
structure of InSb, and noting that the Fermi en- 
ergy on the m side (0.05 - 0.10 ev) is considerably 
greater than that (0.02 ev) on the p side, it is, 
therefore, impossible (in our junctions) for a 
carrier with entirely transverse energy to tunnel 
from the Fermi surface in the conduction band to 
the light-hole band. Thus, the observed oscilla- 
tions must be due to electron-heavy hole tunnel- 
ling, in which case the selection rules are re- 
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laxed and the dependence of the amplitude of the 
oscillations on H will not be simple. 

Some remarks can be made on the absence of a 
second component in the oscillations contributed 
by Landau levels in the valence band. Because 
of instrumental limitations, fields of 15 kgauss 
or more are required for the observation of the 
electron Landau levels. Hence, those for the 
heavy holes would require considerably higher 
fields. Another factor tending to obscure the 
Landau levels on the p side will be the greater 
scattering frequency due to the relatively greater 
impurity concentration. On the other hand, it 
might be possible, in principle, to detect oscilla- 
tions in the tunnel current of electrons from well 
below the Fermi level in the conduction band to 
the Fermi surface in the light-hole band. How- 
ever, these have not been observed, the probable 
reason being that at the fields at which oscilla- 
tions can just be detected (~15 kgauss), the maxi- 
mum energy of the light-hole band is already 
lower than the Fermi energy in the valence band. 

In the case of E ||H, the infinite degeneracy of 
Landau states is responsible for the sharp peaks 
in the density of states. In the E1H geometry, 
this degeneracy is broken since the energy of a 


Landau level now depends upon the position of the 
electron in the plane normal to H. Thus the den- 
sity of states will vary with H in the perpendicular 
geometry, so accounting for the absence of oscilla- 
tions. 

The fact that, to within experimental error, the 
Landau spacings at the Fermi level in the con- 
duction band determine the current oscillations 
at relatively high biases implies that the tunnel 
current is still dominated by transitions to or 
from states close to the Fermi level on the n 
side. 

The authors are grateful to J. K. Galt for pro- 
viding the magnet facility, and to W. L. Feld- 
mann, J. F. Gilbert, and H. W. Dail for assist- 
ance with the experiments. 
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LOW-TEMPERATURE TRANSPORT IN “SPLIT ~-GERMANIUM” 


Seymour H. Koenig and John J. Hall 
International Business Machines Watson Research Laboratory, Columbia University, New York, New York 
(Received November 28, 1960) 


By measuring conductivity and Hall mobility of 
p-Ge at low temperature under large strains 
(>10°*), we have been able to determine the sign, 
the ratio, and the approximate magnitude of the 
two deformation potentials that describe the 
change, under shear, of the valence band. In 
addition we have been able to verify that the ve- 
locity dependence of the hole-acceptor recombi- 
nation cross section is independent of the ac- 
ceptor ground-state energy over a limited range, 
in agreement with the prediction of the trap 
mechanism proposed by Lax.' 

When germanium is subjected to a shear stress 
the valence band, normally 4-fold degenerate 
(including spin) at the Brillouin zone center, 
splits into two doubly degenerate bands separated 
by an energy 5&.?~* In a neighborhood within an 
energy << 56 of each (new) band edge, the sur- 
faces of constant energy are ellipsoids of revo- 
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lution about the stress direction, and the effect- 
ive masses are functions of the geometry of the 
strain, and not its magnitude.*»*,° 

For a uniaxial stress y in either the (100) or 
{111) directions, the splitting is given by 


685 99 = 244 99 = 254 (1 +A) 1 Ox! =4S,,(1 +X) ID xi/3 
=21Dl€, 99 = 2-44x10" 7 1x51 ergs, (1a) 
66554774411 =V8S,, \dx|/3=2S,,|D,'x!/3 
=VB1d1€ 49 1S q4MS1 1 +2549 +844) 
= 0.838 x10°'*|yd| ergs, (1b) 


and b= -2D,,/3; d=-2D,'/V3. Here b is the de- 
formation potential constant (for hole energies) 
appropriate to a strain of tetragonal symmetry 
(with respect to the conventional cubic basis) as 
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defined by Pikus and Bir® (PB); d is the similar 
constant for a deformation of rhombohedral 
symmetry. D, and D,’ are similar constants 
defined (for electron energies) by Kleiner and 
Roth* (KR). The definitions used are those in 

Eq. (2) and Table I of KR, rather than their ver- 
bal description of D,, and D,,’, which is apparently 
inconsistent with these. S,,, S,,, S,, are the 
elastic compliance constants, and A the Poisson 
ratio. €,9, and €,,, are the respective longitudinal 
components of strain that result from the stress 
y.° The relation to the notation of other authors 
is given in Table I of Price and Kao.’ 

Sufficiently large stresses may be applied to 
germanium so that at low temperatures ($10°K), 
5&>>kT. In this limit only one band is appre- 
ciably populated; the Hall mobility then becomes 
independent of strain. Figure 1 shows the strain 
dependence of the Hall mobility for two samples 
of 20 ohm cm p-type Ge cut from the same crys- 
tal, one (111) oriented, the other (100). The 
abscissa scale was chosen so that if b=d, it 
would be proportional to 56 for both samples. 

The saturation of the mobility is clearly indi- 
cated. From these data, a ratio |d/b| =2.1+15% 
or |D,,’/D,,|=3.6 is obtained. This result is a 
factor of ~2 less than the estimates made by PB* 
and KR* from piezoresistance data. It is difficult 
to obtain any precise value for d from the present 
data. However, if one assumes that when the 
mobility is at 90% saturation, 56 ~3kT, then 

d=1 ev. 

According as the sign of b and d, the relevant 
effective-mass ellipsoid is either oblate or pro- 


late.* To determine the signs we have crudely 
measured the transverse conductivity under sa- 
turation strain for both samples by making a two- 
terminal measurement of the resistance between 
the Hall contacts. For the (100) case, pu y/Hy> 1 
corresponding® to b<0 or D,>0, ignoring the 
(unlikely) possibility of an extremely large scat- 
tering anisotropy. Similarly, for the (111) case, 
Hy /uy <1 corresponding to d<0 or D,,’>0. 

These results agree with the signs of the esti- 
mates made by KR‘ and those which may be de- 
duced from PB.° It is interesting to note, how- 
ever, that the work of Hensel and Feher’® yields 
the opposite signs for the corresponding valence 
band deformation potentials in silicon.* (Actually 
only the signs of the products dD and bB are here 
determined experimentally, where D and B are 
the conventional matrix elements of the spin- 
orbit interaction; but the signs of D and B are 
presumed known.®*») 

If 5& can be made sufficiently large, the ac- 
ceptor ground state will be formed from Bloch 
states near the edge of the upper of the two split 
bands. The effective-mass approximation as ex- 
tended to the prolate case by Keyes,° using the 
masses calculated by Pikus and Bir,° and as- 
suming that the dielectric constant remains iso- 
tropic, gives for the saturation activation energy 
€ the value 4.25x107° ev for a (100) stress. 
Figure 2 shows a plot of reciprocal Hall constant 
vs strain for the (100) orientation at T=6°K. The 
maximum stress on the (100) sample, which cor- 
responds to a strain of 0.7x10-?, was limited by 
the available apparatus. (This corresponds to 
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FIG. 1. Hall mobility > 
(normalized to zero , 
strain value) vs strain for 
a (100) and (111) oriented 
sample of p-Ge. The ab- 
scissa scale is in units of 
the strain-dependent split- 
ting of the valence band at . 
the zone center divided by 
the appropriate deforma- 
tion potential [see Eq. (1)]. 
The strains at saturation 1.0 
are >107%, 
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Reciprocal Hall constant vs strain. The 
variation is due to a change in impurity activation en- 
ergy which, in terms of the effective-mass approxi- 
mation, is caused by a large reduction in at least one 
component of the mass tensor. 





5§ ~7x107* ev.) It would appear from the data 
that the condition € << 5& has almost been reached. 
The activation energy measured at this strain is 
6.8x10°* ev. Whether there is indeed a discrep- 
ancy with the predictions of the effective-mass 
approximation remains to be determined. Fora 
(111) stress, the reciprocal mass tensor is 
oblate; the calculated ¢€ in this case is 5.0x10°° 


This “new” semiconductor, “split p-germa- 
nium,” has a unique application in hot electron 
problems for values of strain and electric field 
such that 66>>kT, where T, is the electron 
“temperature.” The nonlinear current-voltage 
relation": should then be independent of strain 
(except for a strain-dependent scale factor in the 
current caused by the varying activation energy) 
until impact ionization of neutral impurities 
makes a significant contribution to the current 
density. The electric field at which this occurs 
has been a matter of conjecture for some time.” 
However, it is now possible, by varying the 
strain, to vary continuously the impurity activa- 
tion energy, without changing any transport 
Figure 3 shows recorder traces of 
































FIG. 3. Variation of current 
with electric field for two values 
of strain sufficiently large so 
that the transport parameters 
are strain independent but the 
activation energy not. The two 
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curves diverge when impact ioni- 
zation of neutral impurities be- 
gins to contribute appreciably to 
the carrier density. (The cur- 
rent scales have been adjusted 
to compensate for the activation 
energy change by making the 
curves overlap in the Ohmic 
range.) 
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current density vs electric field for two large 
values of strain; the current sensitivity has been 
adjusted to compensate for the activation energy 
change (~20%). There is overlap of the two 
curves for a considerable range of electric field 
where the nonlinearity in conductivity is due to 
the velocity dependence of the recombination 
coefficient o.’* It is seen that o varies by a 
factor of 3.5 over this range. The (velocity- 
dependent) recombination cross section is ex- 
pected to be dependent, not on the acceptor 
ground-state energy, but rather on a continuum 
of states ~kT above the valence band edge,’ 
which depends in this case on the parameters of 
the upper band. This expectation is here con- 
firmed experimentally, and is additional strong 
confirmation that the cascade recombination 
mechanism of Lax is operative. 

We are indebted to Rodney Brown for computa- 
tional assistance and to Peter J. Price for dis- 
cussions. 
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POLARIZATION OF THE CONDUCTION ELECTRONS IN THE FERROMAGNETIC METALS* 


A. J. F. Boyle, D. St. P. Bunbury, and C. Edwards 
The Physical Laboratories, University of Manchester, Manchester, England 
(Received November 15, 1960) 


It has recently been proposed’ that the polar - 
ization of the 4s conduction electrons in Fe, Co, 
and Ni is antiparallel to that of the electrons in 
the unfilled 3d shell; this is the reverse of the 
usual assumption. There are few ways in which 
this polarization is manifested. One, however, 
is in the effective magnetic field at a nucleus, 
which acts through the contact interaction with 
the 4s electrons. The magnitude of this field has 
been given by Marshall’ as 


H = (81/3)é || (0) |, Unb, (1) 


where p is the Bohr magneton, ~,1(0)1 4° is 
the average probability density of a 4s conduc- 
tion electron evaluated at the nucleus [| (0)! a 
being the free atom value], m is the number of 
conduction electrons per atom, and p is their 
polarization. The sign of the field is positive, 
i.e., parallel to the direction of magnetization, 
if the polarization of the 4s electrons is parallel 
to the 3d polarization. 


Hanna et al.*»* have measured the magnitude 
and sign of the field at the Fe nucleus in iron by 
observing the Zeeman splitting of the 14-kev 
transition in Fe*’ and obtain the value -3 x10° koe. 
In this case, however, there are other contri- 
butions to the field beside H,; mainly those due 
to the electrons of the same atom, i.e., the 
polarized 3d electrons. These other contribu- 
tions almost certainly outweigh the effect of the 
field H,, alone (H,. is probably of the order of 50 
to 100 koe), and it is therefore difficult to draw 
any conclusion concerning its sign. 

In an attempt to measure H,, directly, dilute 
(1%) solid solutions of Sn in the ferromagnetic 
metals (Fe, Co, and Ni) have been prepared. 
Since Sn is basically diamagnetic, we expect 
that there will be no contribution to the field at 
the Sn nucleus from its own inner electrons and 
that the field will be given simply by 


H , =(41/3)M +H (2) 
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where |¥(0)| 4” in (1) is the appropriate value for 
a Sn atom. Since H, is due to an over-all polar- 
ization of the conduction electrons, the other 
terms should remain those appropriate to the 
solvent atoms.°® 

The Zeeman splitting of the 24-kev transition 
in Sn"*® was observed® using the Méssbauer 
effect.7**> The transmission through the absorber 
(the alloy) of the radiation emitted by a moving 
source was measured as a function of the velocity 
of the source. The velocity spectrometer has 
been described in detail elsewhere.® The tem- 
perature of both the source and absorber were 
maintained around 100°K. Figures 1, 2, and 3 
show the transmission spectra obtained for Fe 
and Ni and Co. Each member of the doublet is 
an unresolved triplet caused by the splitting of 
the excited 3" state, while the doublet separa- 
tion is due mainly to the splitting of the ground 
3* state. The field strengths H. listed in Table I 
were derived using the known values of the mag- 
netic moments. (The excited state moment has 
been measured as 0.83+0.03 nuclear magneton.) 

The shifts listed in the table refer to the dis- 
placement of the spectra towards negative vel- 
ocity and are almost entirely due to the chemical 


velocity mm/sec 


Table I. Results derived from the velocity spectra 
of the absorption of 24-kev y rays from Sn'!* in ab- 
sorbers containing 1% of Sn''* dissolved in Fe, Co, 
and Ni. The shift is the displacement of the spectra 
towards negative velocity due to chemical effects, H, 
is the field at the Sn nucleus, and (47/3)M is the usual 
Lorentz field. 








Shift H, 
Solvent (mm/sec) (koe) (47/3)M H,-(40/3)M 
Fe pe | -81+4 7.5 -88 
Co a -20.5+1.5 6.1 -26.5+1.5 
Ni 1.1 +18.52+1 2.0 +16 .521 





shift.*° In this case, where the source was 
metallic tin, the magnitude of the shift should be 
proportional to (&, - g,)1¥(0)| gy"; Esp and é, 
refer to metallic tin and the solvent metal, 
respectively. 

The sign of H, was determined by placing each 
absorber in a transverse magnetic field of about 
7 koe; the magnetization direction is now parallel 
to the external field. A positive field was ob- 
served only in the case of Ni. The final values 
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FIG. 3. The absorption 
spectrum obtained with 1% 
Sn'!9 in Co. 
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of [H, - (41/3)M] are listed in the table, where 
the values of M have been corrected for the 
quenching of the atomic moments by the added 
valence electrons of Sn. This is most significant 
in Ni, and, since it will occur predominantly 
amongst the nearest neighbors, will undoubtedly 
result in a lowered value of the observed field. 


The results are obviously inconsistent with the 
existence of the single field H,; firstly because 
of the different signs and secondly because of 
the relative magnitudes of the fields. Since p is 
very closely proportional to M,? the expression 
(1) for H,, can be written simply, H,.=const &.M. 
Further, we can conclude from the chemical 
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shifts that £, is practically constant for Fe, Co, 
and Ni, and therefore H,. will be roughly pro- 
portional to M. 

There is another mechanism" by which an 
effective field might be produced at the Sn nucleus 
in this situation. The wave function for the 4s 
electrons of the Sn atom will overlap with those 
of the 3d electrons of the surrounding solvent 
atoms, and the polarization of the latter will re- 
sult in a change in the relative spatial distribu- 
tion of the 4s electrons in respect of their spin 
orientation. In the region of the Sn nucleus, 
electrons with spin antiparallel to the 3d polar- 
ization will predominate, producing a negative 
effective field H,. 

Evidence for the existence of such a field is 
provided by the results for Co (Fig. 3) which 
indicate the presence of a second field of -50 koe 
with the same chemical shift. X-ray analysis of 
the sample showed that both cubic and hexagonal 
structures were present in the rough proportion 
40:60. Since neither £, nor p depends on the 
structure, the two values of the field cannot be 
associated with H.; however, since the distri- 
bution of the 3d orbitals will probably differ for 
the two structures,’ Hp would also differ. 

We have no estimate of H pb» but its magnitude 
should depend not only on M but also on the mean 
radius of the 3d shell and thus will decrease 
more rapidly from Fe to Ni than does the field 
H_,. which depends only on M. In the absence of 
any further contributions, combination of such a 
field Hp with a positive value of H,. would there- 
fore provide a qualitative interpretation of the 
present results. 

It is interesting to compare the present results 
for the field at a Sn nucleus in Fe with those of 
Samoilov et al.’* From measurements of the 
nuclear polarization these authors obtain values 
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of 250 koe and 280 koe, respectively, for the 
field at In*“* and Sb”? dissolved in Fe. 

We are indebted to Dr. C. Johnson for drawing 
our attention to the results of Samoilov et al. 
and we are grateful to many people for valuable 
discussion, particularly Dr. Lomer, Dr. W. 
Marshall, Dr. S. F. Edwards, and Dr. H. E. 
Hall. Much helpful advice in the preparation of 
the alloys has been given by Dr. J. Stubbles. 
We would also like to thank Dr. J. Zussman for 
performing the x-ray analysis. 
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STIMULATED INFRARED EMISSION FROM TRIVALENT URANIUM 


P. P. Sorokin and M. J. Stevenson 
Research Center, International Business Machines Corporation, Yorktown, New York 
(Received November 28, 1960) 


Schawlow and Townes?’ have proposed the ex- 
tension of maser techniques to generate coherent 
light in the visible and infrared regions of the 
electromagnetic spectrum. Recently Maiman? 
and Collins et al.* obtained stimulated emission 
of radiation in chromium-doped Al,O, in the R, 
transition from the *E metastable state to the 
ground state. The oscillation was achieved by 
inverting the population of the metastable state 
relative to the ground state by a sufficient amount 
to overcome the losses in the optical cavity. 
Since this required approximately 10** - 10** ions 
in the metastable state, very high optical pump- 
ing power was necessary to achieve inversion. 

This Letter describes the characteristics of 
stimulated emission from trivalent uranium ions 
substituted for divalent calcium ions in calcium 
fluoride. The crystals were grown by a technique 
similar to that described by Galkin and Feofilov.* 
In contrast to the case of the ruby optical maser, 
the oscillation takes place in a transition from a 
metastable state to a level approximately 515 
cm! above the ground state. At low temperatures 
this state is depopulated by at least a factor of 
10*° relative to the ground state. The general 
condition for stimulated emission is that the num- 
ber of excited ions in the metastable state must 
exceed the number of ions in the lower state by a 
certain amount An sufficient to overcome the 
cavity losses. In trivalent uranium this condition 
will be achieved at low temperatures when the 
population of the metastable state is approxi- 
mately equal to An since the population of the 
lower state is negligible. The value An may be 
generally much less than the number of ions in 
the ground state. It was indeed observed that the 
pumping power necessary to achieve stimulated 
emission was lower by about a factor of 500 com- 
pared to ruby. 

Figure 1 shows the pertinent energy levels of 
the trivalent uranium ion when substituted for 
calcium in the CaF, lattice.‘ The ion has several 
strong absorption bands in the visible and near 
infrared wavelength regions. The strongest ab- 
sorption bands are in the green and blue, and as 
a result the crystals have a pale to deep red 
color and thus show a striking resemblance to 
ruby. Very strong fluorescence is observed in 
the infrared region between 2.0 and 2.6 microns 
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FIG. 1. Energy level diagram of trivalent uranium 
in calcium fluoride based on the paper by Galkin and 
Feofilov. Broad-band pumping light applied in the blue 
and green visible spectrum causes transitions to ex- 
cited bands, one of which is depicted in the figure. 
These pumping transitions are followed by rapid, non- 
radiative transitions to the two metastable levels shown. 


when exciting radiation is applied in the range of 
the absorption bands. The room temperature 
fluorescence spectrum consists of four peaks. 
The two shorter wavelength peaks coincide with 
the absorption peaks and therefore are resonance 
lines. The two longer wavelength peaks repre- 
sent emissions to a state located approximately 
515 cm~* above the ground state. The longest 
wavelength emission occurs from the lower 
metastable state and is the transition in which 
the maser oscillations are observed. 

A single crystal of calcium fluoride with 0.05 
mole percent uranium was cut into a cylindrical 
shape 3/8 in. in diameter and 13 in. long. The 
ends were polished flat to one-twentieth of a 
sodium wavelength and parallel to fifteen seconds 
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of arc. Silver films were evaporated on both 
ends, with one end left completely opaque and 
the other end left with a 1.5% transmission. A 
protective SiO coating was applied over the silver 
to prevent oxidation. The sample was placed in 
a specially constructed liquid helium Dewar and 
was cooled without being in actual contact with 
the liquid helium. A pulsed xenon discharge 
lamp excited the vertically positioned sample 
through a side port in the Dewar. The infrared 
radiation generated in the sample was coupled 
out through the 1.5% transmitting end face and 
was detected by a lead sulfide detector through a 
top port. Either a 2.5-micron interference filter 
or a Perkin-Elmer monochromator was used to 
isolate the infrared radiation from the lamp 
background. The response time of the detector 
was approximately 40 microseconds for decaying 
radiation. 

Figure 2 shows typical oscilloscope traces of 
the detector response to the intense monochro- 
matic stimulated emission when the pump power 
is set above the threshold value. The signal 
during the first 450 microseconds of Fig. 2(a) is 
due primarily to the scattered light from the 


(Qa) 





(b) 





Fig. 2. (a) Detector response above threshold show- 
ing the onset of oscillation after 450 usec. A filter with 
a bandwidth of 0.1 micron and centered at 2.5 microns 
is inserted before the detector. The sweep speed is 
200 usec/division. (b) Detector response with the con- 
ditions as in 2(a) except that the oscilloscope gain is 
reduced by a factor of 500 and the sweep speed is 100 
usec/division. 
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xenon flash lamp plus a small but undetermined 
amount of spontaneous radiation. Oscillation 
begins where the rapidly rising, almost vertical, 
trace is observed and always occurs when the 
pumping power reaches the same critical level. 
Fig. 2(b) shows the same detector output with 
reduced gain and expanded time scale. The ratio 
between the stimulated emission and the back- 
ground light was measured to be in excess of 
1000, which means that the spontaneous emission 
is less than the stimulated emission by a factor 
of several thousand. This ratio is larger than 
that observed for the ruby maser, which con- 
forms to the expectation of smaller spontaneous 
emission in the infrared as compared with the 
visible spectral region. The spectral width of 
the maser pulse was substantially less than the 
resolving power of the monochromator. The 
beam was confined to an angle of approximately 


(a) 


b) 





(Cc) \ 
—— 


FIG. 3. Relaxation oscillations at varying pump 
power levels. The time scale is 100 usec/division. 
Trace (a) shows the maser output with pump power 
about 1.5% above threshold. Traces (b) and (c) are 
taken at powers barely above threshold. 
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0.01 radian. 

An interesting feature of the maser output 
shown in Fig. 2(b) is the presence of relaxation 
oscillations, previously also observed in less 
detail by Collins et al. in ruby. The number of 
observed spikes decreases rapidly as the power 
input into the exciting xenon lamp is lowered 
towards the threshold value until at barely above 
threshold only a single spike oscillation occurs. 
Figure 3 illustrates this phenomenon. The 
single spike takes place at exactly the peak of 
the pulse output of the xenon lamp, while the 
multiple spikes begin earlier because the ex- 
citing power reaches the threshold value at an 
earlier point. The spacing between the oscilla- 
tions increases as the xenon lamp input power is 
brought closer to threshold. These features seem 
consistent with the explanation offered by Collins 
et al. for the relaxation oscillations in the ruby 
optical maser, perhaps complicated by contribu- 
tions from several frequency modes within the 
inhomogeneously broadened line as suggested by 


G. J. Lasher of this laboratory. 

The favorable energy level structure of triva- 
lent uranium as compared to ruby may result in 
a strong cw infrared source for research pur- 
poses in the near future. 

We would like to thank J. Lankard for assist-.- 
ance with the experiments, W. V. Smith and 
A. H. Nethercot for their continued interest and 
encouragement, and other members of the labo- 
ratory for helpful discussions. Special thanks 
are due Walter Hargreaves of Optovac, Inc. for 
cooperating with us in growing the crystals. 
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MOBILITY OF HYDROGEN IONS* 


Robert N. Varney 
Department of Physics, Washington University, St. Louis, Missouri 
(Received November 14, 1960) 


Experimental measurements of the drift veloc- 
ity or mobility of hydrogen ions in the parent gas 
indicate only a single type of ion, which has con- 
sistently been assumed to be H,*.'~* Rarely, a 
second ion has been reported,*** but it has always 
been difficult to observe so that a cloud of suspi- 
cion as to its reality has prevailed. Identifica- 
tion of the ion as H® seems to be justifiable. The 
theory of ionic mobilities has been applied by 
Mason and Vanderslice® to hydrogen, and it sup- 
ports the identification of the ions described 
above. It predicts also that H,* should be found 
with a mobility considerably higher than that of 
either H* or H,* or of any observed values re- 
ported in references 1-5. 

Meanwhile, Stevenson’ has published some 
properties of the H,* ion which are listed here: 

(1) H,* is formed by the reaction H,* + H, = H,* 
+H. 

(2) The cross section for this reaction is enor- 
mous, approximately 110°" cm. 

(3) The binding energy of H,” is 97 kcal or 
4.18 ev (private communication). 


(4) H,* ions passing through H, gas probably 
have an enormous charge-exchange cross section, 
of the same approximate size, 1x10~* cm’, as 
for the formation of H,". 

(5) The charge-exchange mechanism appears 
to be an exchange of a proton instead of the more 
usual exchange of an electron envisaged by the 
term “charge exchange.” 

If these findings are correct, and there is 
every reason to believe they are, then two 
changes must be made in the theory and inter- 
pretation of mobilities of ions in hydrogen. 

(1) The H,* ion is never observed in mobility 
experiments nor will it ever be observed. 

(2) The .” ion has a lower mobility than pre- 
dicted by Mason and Vanderslice, presumably 
having the value commonly observed. 

The following justifications are offered: The 
chemical equlibrium constant for the formation 
reaction for H,* given by Stevenson is so large 
that there is no possibility of shifting the equili- 
brium toward the H,* side of the reaction within 
the bounds of swarm-type experiments, for ex- 
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ample, by very low pressures. At very high 
temperatures or very high values of E/p,, at 
which equilibrium would be shifted toward the 
left side, the H,* ion would still not be found as 
the dissociation of H,* would then prevail. It 
may be noted that the binding energy of H,* at 
2.65 ev is appreciably less than that of H,* quoted 
as 4.18 ev by Stevenson. 

The structure of H,* seems necessarily to be 
H-H-H. Since H, is simply H=H, the exchange 
of the proton is not surprising. The binding en- 
ergy of H, is very nearly the same as for H,", 
4.48 and 4.18 ev, respectively, so that the proton 
exchange is virtually a resonance phenomenon. 
The H,” ion passing through H, gas is, in effect, 
the “normal” ion of the parent gas consisting of 
the parent molecule with an attached proton in- 
stead of the parent molecule with a detached 
electron. The size of the charge-exchange cross 


section should drastically lower the ionic mobil- 
ity. 

The writer is much indebted to Dr. Stevenson 
for pointing out the experimental evidence he had 
obtained and for extensive discussions in which 
the analysis above was formulated. 
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DISSOCIATION OF H,* AND He~ BY ELECTRIC FIELDS 


A. C. Riviere and D. R. Sweetman 
United Kingdom Atomic Energy Authority, Aldermaston, Berkshire, England 
(Received November 28, 1960) 


Processes leading to the dissociation of mole- 
cular ions have assumed some importance re- 
cently in connection with injection into thermo- 
nuclear devices. In particular, the dissociation 
of molecular ions by electric fields (or the equi- 
valent ¥xB field), provides an attractive method 
of injection of high-energy particles which avoids 
the use of dense low-temperature plasmas or 
carbon arcs. 

Recent theoretical studies by Hiskes et al.'~* 
have shown that the upper vibrational states of 
molecular ions may dissociate in electric fields 
of a few hundred kilovolts/cm. In this paper we 
describe experiments in which this dissociation 
has been observed and also an experiment in 
which the He ion was dissociated by an electric 
field. 

The experimental arrangement was basically 
that described by Sweetman‘ for measuring the 
H,* dissociation cross section in hydrogen. For 
the present measurements the gas target cham- 
ber was replaced by an electrode system so de- 
signed that the ion beam would pass through an 
electric field of known strength. Figure 1 shows 
the details of the high-field gap and the field dis- 
tribution as seen by a particle on the axis and on 
the periphery of the beam. The ions passed 
through the gap and the high-voltage electrode, 


560 


HIGH 
VOLTAGE 
SUPPLY 


ro 
( 








(ON_BEAM _ p 


nae an ae 


A\ ims 
ar 



































4 P 
0.24 MM DIA. 106mm 
we .. ELECTRIC FIELD 
1.oF : = IN IDEAL CASE 
fey) OFF. 
sas ON AXIS 
0.6} 
O.4F 
MODIFIED FIELD 
DUE TO PRESENCE 
7 OF APERTURES 
re) L 1 
-.OMM  _—O.5MM re) +0.5MM +1.OMM 


FIG. 1. Details of the high-field gap and electric 
field distribution. 
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and then returned to ground potential through a 
region in which the field was everywhere less 
than 2% of that in the gap. The ion beam was 
magnetically analyzed and the individual com- 
ponents were detected by CsI scintillation count- 
ers. Coincidences between H° and H* particles 
and also between two H* particles were recorded 
so that both modes of dissociation could be 
studied under low-background conditions. 

In the first experiments, H,* ions emerging 
from an rf ion source were accelerated to 2 Mev 
in a Van de Graaff electrostatic accelerator be- 
fore passing into the apparatus. The fraction of 
the beam dissociated by electric fields up to 
510° volts/cm is shown in Figure 2. The num- 
ber dissociated increases steadily for fields 
above 0.70x10° volts/cm. The rate of ionizing 
dissociation (to H* + H*) remains constant over 
this range. This is to be expected since the field 
strength required to ionize hydrogen atoms is of 
the order of 10° volts/cm,! whereas the maxi- 
mum field strength achieved here was 510° 
volts/cm. This result also removes any doubt 
about possible dissociation due to the deflection 
of the ion beam on to stops by the applied elec- 
tric field. No such effect is observed. 

In order to increase the population of the upper- 
most vibrational levels, H,* ions were obtained 
from H,* ions which had been accelerated to 2.0 
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ELECTRIC FIELD IN KILOVOLTS/CM 


FIG. 2. Dissociation of H,” ions as a function of 


electric field. Curves (1): H,” ions direct from source. 


Curves (2): H,* ions obtained from the breakup of H;* 
ions. 


Mev and then dissociated by collision with gas 
molecules. It can be seen that the H,* ions 
formed by the breakup of H,* ions were dissocia- 
ted by the field to a much greater extent than 
those obtained directly from the ion source. This 
mode of formation evidently results in a con- . 
siderably greater population of the upper vibra- 
tional states. 

The calculations of Hiskes* show that the upper- 
most vibrational levels (v=17, 18) with zero ro- 
tational energy should be completely dissociated 
by fields of 6x10° and 210° volts/cm, respect- 
ively. A steep rise in the fraction dissociated 
would be expected at 2 x10° volts/cm followed by. 
a plateau. In fact the dissociation fraction is 
seen to rise almost linearly from fields of about 
0.710° volts/em. According to Hiskes,’ the 
presence of rotational energy probably has the 
effect of reducing the field required for dissocia- 
tion, and this may account for the lack of distinct 
steps due to the discrete vibrational states. 

The dissociation of HD*, D,* and other mole- 
cular ions is of great interest for thermonuclear 
experiments and will be investigated in future 
work. 

The effect of an electric field was also tried 
on He and H™ ions at energies of 1.0 Mev and 
0.33 Mev, respectively. The results are shown 
in Fig. 3. The He ion has an estimated binding 
energy of 0.075 ev,° and was completely disso- 
ciated by a field of 4.5x10° volts/cm. The H™ 
ion has a binding energy of 0.75 ev,° and no effect 
was observed up to the maximum field of 5.4 10° 
volts/cm. Holgien and Midtdal® have shown that 
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FIG. 3. Dissociation of He” and H™ as a function of 
electric field. 
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the probable configuration of He™ is (1s 2s 2p)*P, 
and that autoionization from this state is not al- 
lowed. In the absence of an electric field the 
lifetime of this state is known to be greater than 
10-* second.” The dissociation of this ion in the 
electric field may be due to breakdown of the 
selection rules for autoionization or due to field 
emission of the electron leaving the He® in the 
(1s 2s) metastable state. Further experimental 
investigations are under way to elucidate this 


point. 
The authors are grateful to Mr. R. E. Brad- 


ford, Mr. J. Coupland, and Mr. R. Webster for 
their help in setting up the apparatus and taking 


the results. 
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SEARCH FOR A SMALL CHARGE CARRIED BY MOLECULES* 


John G. King 
Department of Physics and Research Laboratory of Electronics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 
(Received November 7, 1960; revised manuscript received November 28, 1960) 


The suggestion by Lyttleton and Bondi’ that 
hydrogen atoms might carry a charge of approx- 
imately 2x10~** of an elementary charge has 
stimulated a number of investigators to under- 
take experiments designed to measure this charge. 
At least two methods have been used: atomic 
beam deflection experiments,”** and the method 
of Piccard and Kessler.* The results published 
thus far?~® are listed in Table I. In this paper 
preliminary results of a Piccard-Kessler exper- 
iment are reported in which hydrogen molecules 
(and also helium atoms) exhibit a charge that is 
approximately 40 times less than that required 





by the hypothesis of Lyttleton and Bondi. 

In our experiment the gas whose charge is to 
be measured is allowed to escape from an elec- 
trically insulated metal container connected to an 
electrometer. While the gas is escaping, a cur- 
rent that can be measured will flow to the con- 
tainer if the gas is charged. The charge on each 
atom can then be computed, although it is essen- 
tial to distinguish charge carried by each of the 
hypothetically charged atoms from charge carried 
by ions, electrons, and charged dust particles. 
By passing the escaping gas through a de-ionizer 
consisting of a coaxial capacitor charged by a 








Table I. Upper limit on the charge carried by various molecules. 
Charge 
Molecule (elementary charge/molecule) Experimental method 
co, < 2x107'% Piccard-Kessler 2 
CsI < 4x1078 Molecular beam 
A < 8x 10770 Piccard-Kessler © 
N, < 12x 107?° Piccard-Kessler © 
CsF <= 2x10-" Molecular beam 4 
KF < 1x107% Molecular beam4 
n < 6x10-" Molecular beam & 








dsee reference 3. 
©See reference 6. 


4See reference 4. 
bgee reference 2. 
°See reference 5. 
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battery within the container, ions and electrons, 
because of their relatively high mobilities, can 
be swept from the gas stream, while atoms 
carrying the hypothetical charge of 107” of an 
elementary charge will be unaffected. Dust 
particles, which range in size from molecules 
to visible grains, may be sufficiently massive to 
pass through the de-ionizer, but their effect is 
likely to be both erratic and cumulative. 

By recording the pressure and the output of 
the electrometer, both as a function of time, on 
a two-channel recorder, it is not only possible 
to measure the total charge carried by the gas, 
but also to see that the potential-versus-time 
curve has the expected shape for the flow of a 
uniformly charged gas. The course to be fol- 
lowed by the experimenter is plain. If the atoms 
seem to be charged, the computed charge per 
atom must be independent of all parameters, 
such as de-ionizer potential and polarity, rate 
of gas efflux, pressure, bias potential of the 
container with respect to ground, and, of course, 
details of apparatus construction and gas hand- 
ling. It is to be expected that the charge per 
atom will be different for different gases. On 
the other hand, if the atoms appear uncharged, 
it is important to make certain that the electrom- 
eter is sensitive during the gas efflux and that 
mechanisms that might cancel a genuine effect 
are eliminated. Such cancellation is most likely 
to occur at the end of the outlet pipe from the 
container. 

Figure 1 is a scale drawing of the apparatus. 
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FIG. 1. Scale drawing of the apparatus. 


The outer can serves both as a vacuum envelope 
and as an electrical shield for the container 

(11.2 liters) which stands on three quartz sup- 
ports and is filled with gas through a quartz pipe. 
This quartz pipe is metalized at its ends so that 
it can be soft-soldered into place. Although it. 
was designed to stand a pressure of 1000 psi, 

the container was usually operated at a pressure 
of 265 psi. The outer can was pumped to 107° 
mm Hg to prevent electrical leakage and to re- 
duce undesirable thermal effects. The central 
electrode of the de-ionizer is maintained at +90 
volts, or zero volts with respect to ground, as 
deter mined by the setting of a switch. The ap- 
proximate contact potential to ground of the con- 
tainer can be determined by a movable vane. The 
gas that is to be investigated is admitted to the 
container through a porous bronze filter with 5-y 
pores, an electrostatic precipitator, and needle 
valves for controlling the flow. A pressure 
transducer generates a signal proportional to the 
pressure in the container which can be recorded 
on one channel of a two-channel Sanborn re- 
corder. The gas escapes from the container 
through an external regulator that maintains con- 
stant pressure on a needle valve adjusted to give 
the desired rate of flow. It was thus possible to 
obtain a nearly linear drop in pressure in the 
container with respect to time. 

Potential changes of the container with respect 
to the outer can are measured with a GL5674 
electrometer whose output is amplified by a 
chopper amplifier and recorded on the other 
channel of the Sanborn recorder. Potential dif- 
ferences of 10~* volt, or currents of approxi- 
mately 10° electrons per second, can be detected. 
The time constant is 60+5 sec which agrees with 
the time constant predicted from independent 
measurements of the electrometer grid resis- 
tance and capacitance, 10°? ohms and 60 uu, 
respectively. The ground return of the electrom- 
eter is connected through a low-impedance voltage 
supply that makes it possible to place the con- 
tainer at different potentials with respect to the 
outer can. When the potential of the container 
was +0.01 volt, moving the vane was found to 
have negligible effect. Frequent calibrations of 
the electrometer were made with signals close 
to the minimum detectable signal introduced at 
the ground end of the electrometer grid resistor. 
No consistent dependence of the electrometer 
sensitivity on the polarity of the calibrated sig- 
nal, the amount of gas in the container, the 
polarity of the de-ionizer, or any other param- 
eter was found. Calibrations were also made 
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while the gas was flowing from the container, 
and the sensitivity was found to be normal. 

The 21 hydrogen and 19 helium runs at 265 psi 
and 140 sec efflux time (selected only because 
these parameters were the same) yielded elec- 
trometer recordings that are superimposed in 
Figs. 2(a) and 2(b) normalized at the start of 
flow. Average curves that were constructed for 
hydrogen and helium are plotted in Fig. 3. The 
shaded rectangle, 0.3 mv high and 140 sec long, 
represents 2.7x10° elementary charges, which is 
the amount of charge that would be observed if 
each of the 5.4x10** molecules carried 5x10~*° 
of an elementary charge. The apparent charge 
per molecule, determined from the average 
curves, is (2.5+1.5)x10~*° elementary charge/ 
molecule for hydrogen, and (-4+2)x10~*° ele- 
mentary charge/molecule for helium. The charge 
on hydrogen appears to be negative; the charge 
on helium, to be positive. 

None of the recordings looks like the computed 
curve, (d) of Fig. 3, which was derived from 
the pressure curve, (a) of Fig. 3. The hydrogen 
curve goes negative and then positive, and the 
helium curve has a negative maximum that is 








TIME(SEC) 





~ TIME(SEC) 





FIG. 2. (a) Superimposed electrometer recordings 
of 21 hydrogen runs. (b) Superimposed electrometer 
recordings of 19 helium runs. 
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too early by 30 sec. If it is assumed that the 
helium effect is spurious and that whatever pro- 
duces it is likewise present in the hydrogen runs, 
the helium curve may be subtracted from the 
hydrogen curve. This procedure yields a curve 
(represented by circles in Fig. 3) that is in re- 
markable agreement with the computed curve. 

If this agreement is taken to be significant, we 
obtain a value of (7+ 2.5) x10~-*° elementary 
charge/molecule for hydrogen, and (0+ 2) x107*° 
elementary charge/molecule for helium. 

The electrometer deflections recorded are 
still not understood; but the nature of these re- 
sults, as well as the results of numerous side 
experiments, eliminates a few of the possible 
causes for the deflection, such as capacitative 
effects, piezoelectricity, and thermoelectricity. 
New and improved apparatus, use of other gases, 
including deuterium, and wider variations of 
parameters will provide the answer. In any 
case, the charge that we found for the hydrogen 
molecule, by interpreting our experiment in the 
obvious way, is far less than that required for 
the hydrogen atom by the hypothesis of Lyttleton 
and Bondi. A molecular-beam experiment is 
under way with cesium, but the difficulties of 
achieving sufficient sensitivity with hydrogen 
(and later atomic hydrogen) are considerable. 

My thanks are due to Professor Jerrold R. 
Zacharias and Professor E. M. Purcell for val- 
uable discussions, and to Mr. Frank O’Brien 
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FIG. 3. (a) Observed pressure (arbitrary units) in 


the container as a function of time as the gas flows 
out. (b) Average voltage curve constructed from the 
data of Fig. 2(a). (c) Average voltage curve con- 
structed from the data of Fig. 2(b). (d) Electrometer 
response to the flow, according to curve (a), of a uni- 
formly charged gas from the container. Circles 
represent the difference between curves (b) and (c). 
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for much assistance with the apparatus. 
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CLUSTER MODEL INTERPRETATION OF THE ISOTOPIC SPIN SELECTION 
RULE IN CERTAIN NUCLEAR REACTIONS 


G. C. Morrison 
Enrico Fermi Institute for Nuclear Studies, University of Chicago, Chicago, Illinois 
(Received November 9, 1960) 


In this note we wish to present results concern- 
ing the operation of the isotopic spin selection 
rule in nuclear reactions which have resulted 
from a study of some lithium-induced reactions. 
Experimental verification of such selection rules 
has previously been obtained from the study of 
(d,a)* and (a,d)*? reactions, from the study of 
inelastic deuteron scattering,® and recently from 
the work of Halbert and Zucker* on the inelastic 
scattering of N** by C’*. In the present work we 
have been concerned with the reactions Li®(Li®, 
d)B*° and Li®(Li’, ¢)B*° with particular reference 
to the reactions leading to the T=1, second ex- 
cited state of B’° at 1.74 Mev. In the first reac- 
tion it was to be expected, since both Li® and the 
deuteron have T=0, that the probability of deu- 
teron emission leading to the T=1 level in B’° 
would be reduced by virtue of the isotopic spin 
selection rule; the second reaction was studied 
for purposes of comparison. 

In Fig. 1 the spectrum of deuterons observed 
at a laboratory angle of 9.5° from the reaction 
Li®(Li®, d)B*° is presented. The positions of the 
levels in B’° are as calculated from the observed 
position of the deuterons leading to the first 1+ 
level in B*° and the variation of its position with 
the angle of observation. As can be seen from 
the observed spectrum, deuterons leading to the 
T=1 level in B® are very much inhibited, and 
over the angular range where one could observe 
this deuteron group, the situation persisted. As 
discussed above, this result can be interpreted 
as confirmation of the isotopic spin selection 
rule. However, it is the purpose of this note to 
suggest that the observed inhibition may be con- 
sidered on a different basis. 
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FIG. 1. The spectrum of deuterons observed in the 
reaction Li®(Li®,d)B"® at a laboratory angle of 9.5° 
and incident energy of 2.1 Mev. The positions of the 
corresponding levels in B® are indicated above. 


As part of the reaction studies using lithium 
ions, evidence has been obtained that both the 
reactions (Li®,d) and (Li’,¢) are favored and 
proceed in part by a stripping mechanism with 
capture of the alpha-particle substructure.® In 
the reaction Li®(Li®,d)B*°, if only the J” of the 
final states in B*° are considered, there is no 
reason why the 0+ level should be reduced except 
possibly by a (2J+1) weighting factor. However, 
we see from Fig. 1 that the deuterons leading to 
the 3+ level in B’° are reduced also relative to 
those leading to the 1+ levels. If now the situa- 
tion is considered from the viewpoint of B*° 
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being formed as a result of the capture of an 
alpha particle by Li®, the following considera- 
tions apply. The formation of B’® in a 1+ level 
can proceed by /=0 capture of the alpha particle, 
while the formation of the 3+ level can only pro- 
ceed by at least /=2 capture and might be ex- 
pected to be reduced as is observed. And now 
when we consider the formation of the 0+(T = 1) 
level of B*® on this basis, there is no way by 
which an alpha particle (J =0+) can be captured 
by Li®(J=1+) to form a 0+ state. 

Up to this point, the discussion has only been 
an alternative way of explaining the inhibition of 
the deuteron group leading to the T=1 level, 
since the usual one concerning isotopic spin con- 
servation leads to the same conclusion. How- 
ever, in the comparison study of the reaction 
Li®(Li’, )B*® essentially no tritons leading to the 
T =1 level were observed. In the reaction, there 
is now no restriction due to isotopic spin con- 
servation which can inhibit the formation of the 
T =1 level in B*®. Likewise there is no considera- 
tion due to the available angular momentum for 
the reaction to proceed. 

The spectrum of tritons observed at 14.5° from 
the reaction Li®(Li’, ¢)B’° is shown in Fig. 2, 
with the calculated positions of the levels in B*° 
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FIG. 2. The spectrum of tritons observed in the 
reaction Li®(Li’, t)B'® at a laboratory angle of 14.5° 
and incident energy of 2.2 Mev. The positions of the 
corresponding levels in B’® are indicated above. 
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shown above. The similarity in the relative in- 
tensities of the tritons and deuterons leading to 
the same low-lying levels in B’° can be noted. 
That the second 1+ level in B’® is not more 
strongly excited may be due in part to the low Q 
value (-0.16 Mev) for the Li®(Li’, t)B’°® reaction 
leading to this level. With such small energy 
available in the center-of-mass system the trans- 
fer of an alpha particle to Li® may be inhibited by 
the Coulomb field, although it is realized that 
this argument also may account for some of the 
reduction observed in the tritons leading to the 
0+ level where the Q value is only 0.25 Mev. A 
clarification of this point awaits detailed calcula- 
tions involving the proposed reaction mechanism 
involving alpha-particle transfer. At this time, 
however, the most likely explanation of the inhi- 
bition of the triton group leading to the T=1 
level appears to be the impossibility of Li® cap- 
turing an alpha particle to form a 0+ state of B*® 
as discussed for the Li®(Li®, d)B*° reaction. 

Of course this description of the reaction is 
just what the T=1 nature of the state implies 
insofar as this state in B’® corresponds to a nu- 
cleus with six neutrons and four protons or vice 
versa, where the two extra neutrons or protons 
are necessarily in a singlet state. Thus if one 
thinks of the T=1 level of B*° at 1.74 Mev as 
having the character of Be® plus a deuteron in 
the singlet T=1 state which is consistent with 
the 2-Mev difference in binding energy of the 
singlet and triplet states of the deuteron, the for- 
mation of this level when Li® captures an alpha 
particle would require a rearrangement of the 
deuteron in the target nucleus Li®. That it is the 
T =1 state of the deuteron which is involved in 
both reactions studied and not a T=1 state 
arising from a Li’ + He® configuration can be seen 
from considerations of the nuclear masses in- 
volved which would predict a T=1 level of this 
nature at much higher excitation energies in B’°. 

It would seem therefore that such heavy-parti- 
cle stripping considerations together with such a 
“cluster” model® of the nucleus as discussed here 
are capable of explaining the observed isotopic 
spin conservation in the Li®(Li®,d)B*° reaction in 
a very simple way. Furthermore one sees that 
in addition to the (Li®,d) reaction, the (Li®, a) 
reaction and the (a,d) and (d,q@) reactions, at 
energies where heavy-particle stripping involving 
the transfer of a deuteron is expected to play an 
important reaction role, also fit into this scheme, 
since the observation of the isotopic spin selec- 
tion rule involves a 0+ initial and final state and 
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the deuteron has spin 1. In addition, the conser- 
vation of isotopic spin observed in the inelastic 
scattering of deuterons, alpha particles, or 
more recently C™ ions from N“™ can be under- 
stood as arising from the improbability of the 
deuteron in the target nucleus undergoing an in- 
ternal spin flip to form its T=T state. 

It would appear then that at energies where 
stripping is expected to be important, one does 
not require the introduction of the formal con- 
cept of isotopic spin conservation to describe the 
inhibition observed in reactions leading to such 
0+, T=1 levels. That isotopic spin selection 
rules are experimentally satisfied at high ener- 
gies has also been pointed out by Wilkinson,’ who 
suggests that the reason for their validity is the 
overlapping of many levels in the compound sys- 
tem which break up before isotopic spin mixing 
occurs. In this connection it would be of interest 
to examine the situation in a reaction leading to 
a higher T=1 level where J #0, to see whether 
the transition is allowed or not. It will also be 
of interest to carry out a study of the reaction 
B’°(Li®, d)N***(T = 1) and B’°(Li’, t)N“**(T = 1) to 


verify further the results obtained here. It can 
be pointed out that a study of the latter reaction in 
conjunction with the reaction B*°(Li’, He*)C'*(g.s.) 
should also provide a sensitive test of the rela- 
tions between the cross sections predicted by the 
isotopic spin selection rules. 

The author is grateful for the aid given by 
M. N. Huberman in carrying out the work, to L. 
Palmer and J. Erwood for dependable operation 
of the Van de Graaff accelerator, and to S. K. 
Allison for valuable discussions. 
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SEARCH FOR RESONANCE IN 7-7 INTERACTION IN 1-N SCATTERING AT 0.96 Bev 


J. G. Rushbrooke and D. Radojici¢é 
Cavendish Laboratory, Cambridge, England 
(Received November 21, 1960) 


Photographs of interactions in a hydrogen 
bubble chamber exposed to the 0.96-Bev negative 
pion beam from the Brookhaven Cosmotron were 
kindly supplied by Professor Steinberger and 
have been measured and analyzed; This Letter 
reports a detailed examination of 110 events of 
the type 


Tt +p-p+n +7° 


for evidence of a resonance in the pion-pion in- 
teraction. 

Other photographs of the same exposure have 
been measured by Alles-Borelli, Bergia, Fer- 
reira, and Waloschek’ at Bologna, and some by 
Pickup, Ayer, and Salant® at Brookhaven. Center- 
of-mass momuntum spectra of the 7 and 7° 
mesons from the above reaction were found*** 
to be in reasonable accord with the extended iso- 
baric nucleon model of Lindenbaum and Stern- 
heimer,* which includes contributions from the 
T’= 4 isobaric state (corresponding to the T = 3 


resonance in the 7-N system at T, = 600 Mev) as 
well as from the T’ =3/2 state. 

Laboratory kinetic energy spectra of protons 
were found to possess a low-energy peak which 
has been interpreted by Bonsignori and Selleri® 
and Derado® as evidence for a strong pion-pion 
interaction. This interpretation is based on an 
estimate of the “pole part” of the pion-nucleon 
scattering cross section by the method of Chew 
and Low.’ Their formula, valid when the square 
of the four-momentum transfer A? to the specta- 
tor nucleon approaches the unphysical limit -?, 
where y is the pion mass, becomes in the case 
of the reaction considered, 


870 fA? /u? 1 


sf. 12 
@A2aw" 2m (A7 +7)" q, 


= w( qw* - u?) 
L 





0 (w), 
11 
(1) 


where f? is the renormalized pion-nucleon coup- 
ling constant, q;, the laboratory momentum of 
the incident pion, and o,,,(w) the total cross sec- 
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tion for scattering of the two pions at total energy 
w (in their barycentric system). Bonsignori and 
Selleri argue that if Ong is large enough the pole 
may be expected to dominate the first part of the 
physical zone, and suggest that formula (1) con- 
stitutes a reasonable approximation to the true 
matrix element in the region A?<0.3M’, where 
M is the nucleon mass. 

Upon integration over w*, formula (1) becomes 


80 _ - A? 1 
OA" 2m (A? +y7P hs 
 nax(“") 
xf Mg (ww? /2ut\(1 -4u /u)au?, 
4u 71 


(2) 


where Wmax(4’) is the upper limit obtained from 
the phase space diagram in the (A?, w”) plane. 
Figure 1 shows experimental spectra of our re- 
sults separately, and combined with those of the 
Bologna group, compared with predicted curves 
of* 80/8T from the statistical and isobar models,*® 
and from Eq. (2) on the assumption that o,, is 
constant. The striking low-T peak confirms that 
interactions of the incoming pion with a virtual 
meson emitted by the nucleon are significant in 
the low-A? region. 

Qualitatively one would expect the observed 
low-energy spectrum to be only partially ex- 
plained by this simple mechanism. Peierls*® has 
recently predicted the inelastic 7” -p branching 
ratios on the basis of a model in which final 
states are dominated by two types of final-state 
interaction: (a) between the two final-state pions 


in a state with ¢,,=1, and (b) between a pion and 
the nucleon through the 33 resonance. According 
to this model, about two-thirds of the ( p2~7°) 
cross section involves mechanism (a). This can 
be regarded as an upper limit to the contribution 
of the simple picture described above. This 
agrees qualitatively with the point of view of 
Selleri,! who suggests that the reaction is initi- 
ated wholly by pion-pion interaction in the cloud, 
with a certain probability of their having no 
further interaction with the nucleon, which in the 
(pa~7°) case is probably greater than 50%. 

It is much more difficult to infer the presence 
of a resonance in o,, from available data. Such 
a resonance in o,, would be expected to produce 
a hump in a spectrum of events against w*/”. 

A strong peak at w*/y? = 22 is evident in the 
spectrum of all our events (solid histogram of 
Fig. 2), in disagreement with the predicted curves 
of the statistical and isobar models. Since iso- 
bar formation is known to occur,'~> it should be 
noted that this would yield a general background 
in w*/u*, above which our peak is still signifi- 
cant. Pickup, Ayer, and Salant* report a strong 
peak in their Q distribution for the (1~7°) system 
at Q=325 Mev (w*/y? =19) after considering 
events with protons going backwards in the center- 
of-mass system. 

It is proposed that valuable information on the 
existence and position of a resonance in Ong Could 
be obtained by resolving the available data into a 
number of separate regions of w?/yu?, and con- 
sidering spectra over A?/y? in each. This has 
been done and Fig. 3 compares experimental re- 
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FIG. 1. Laboratory 
kinetic energy spectrum 
of the protons from the 
reaction 1-+p—a~+n°+p 
at incident 1~ kinetic en- 
ergy 0.96 Bev. The histo- 
grams consist of Cam- 
bridge results alone (bro- 
ken lines) and Cambridge 
and Bologna results com- 
bined (solid lines). The 
solid curve is obtained 
from the formula of Chew 
and Low [Eq. (2)], and the 
dashed curves give the re- 
sults of the isobar model 
and the statistical model. 
The theoretical curves 
and the solid histogram 
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FIG. 2. Spectra of Cambridge events against w?/?, 
where w? is the square of the total energy of (~1°) in 
their barycentric system. The solid curves, derived 
from (a) statistical model and (b) isobar model, are 
normalized to the same area as the solid histogram 
which includes all events. Dashed curves were obtained 
by integrating the Chew-Low Eq. (1) over A? up to A? 
=154", assuming (c) 0, to be constant and (d) the 
Frazer-Fulco (reference 13) expression for 0,,(w), 
having w,*=8?, T=0.4; they are normalized to the 
area of the dashed histogram of events with A?<15y?. 


sults with the predictions” of the statistical mod- 
el and of Chew and Low. The phase-space den- 
sity for the Naz system is 


dp = 5( p,? + u”)5(q,” + u?)5(k,” +M”) 
x O( po +92 +2 - DP, - k,)dk,dp,dq,, 


where k, and k, are four-momenta of initial and 
final nucleons, and q,, p2, q2 of the initial pion 
and two final pions, respectively. On integrating 
over superfluous variables, this becomes 


1 


dp = ——_——-(1 - 4? /w?)*dw"da?, (3) 
a,» 
and therefore in the statistical model, 
a Wax (4) 
30; =const J *  (1-4?/w*)7dw*. (4) 
“min 


The strong low-A? peaking of the curves in each 
region obtained from Eq. (1) when integrated be- 
tween the corresponding w* limits is character- 
istic of the A?/(A? +”)? term which has a sharp 
maximum at A?=y?. This peaking is more 


marked than in the total spectrum against T 
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FIG. 3. Spectrum of Cambridge events divided into 
strips of w*/u?; the abscissae are expressed in units of 
A*f?, A solid curve was obtained by integrating the 
Chew-Low formula [Eq. (1)] between the limits of w? 
shown, after inserting the Frazer-Fulco expression 
(reference 13) for o,,(w), with w,?=gu?, [=0.4. A dot- 
dashed curve is the same except that o,, is assumed 
constant. A dashed curve gives the result of the sta- 
tistical theory. Curves of each type are normalized 
to make the area under the total histogram of all events 
equal to the area under the curve obtained by integra- 
ting the corresponding expression over all values of w?, 
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(Fig. 1) and should be easier to recognize in 
practice. Recalling the possible peak at w*/? 
= 22 referred to above, there is no apparent 
resonance effect in the fourth w* “strip” 
(20 <w*/u?<25), the fairly even distribution 
in A? being suggestive only of a two-body break- 
up into nucleon and “di-pion,” to which the special 
assumptions of the Chew and Low theory are not 
appropriate. The low-A? peak in the lowest 
“strip” of w? (4 < w*/u?<10) is perhaps signifi- 
cant, despite the meager statistics, and suggests 
the possibility of a resonance in or near this re- 
gion of w*. Correspondingly a set of curves are 
shown incorporating into Eq. (1) the Frazer and 
Fulco’® ¢,,,=1 resonance expression, in this case 
with w,? =8u?, r=0.4,'*45 

The effect of such a resonance on the spectrum 
of events against w?/y? can be estimated by in- 
tegrating Eq. (1) over A? up to A? =15u?=0.3M? 
from the lower limit given by the allowed region 
in phase space, i.e., 


80 | aa w* 2 /,.2\u/2 
= -4 ( 
dw" 2n45 7 o,0n 2u? alta all 
2 2 
4 ax (o") A? 


x dA’, (5) 


2 rr rer 
4 in (A* +") 


and inserting the Frazer-Fulco Ong ()- The re- 
sulting curve (d) can now be compared with the 
spectrum of events (broken histogram of Fig. 2) 
corresponding to a region in which the pion-pion 
interaction is known to be important from Fig. 1. 
A curve (c) obtained from Eq. (5) assuming Ong 
to be a constant is also shown. The experimental 
evidence is inconclusive when looked at from 
this point of view. 

It is difficult to see how a single Frazer-Fulco 
resonance of narrow width (i.e., [ =0.4) without 
a significant nonresonant contribution could fit 
all the above spectra, as the bulk of events would 
tend to be concentrated into a region of w? near 
w,*, as may be seen in Fig. 2, curve (d) and from 
the solid curves of Fig. 3. Final-state inter- 
actions would tend to mask this effect, however. 
A quantitative test of this point could only be 
made by performing an experiment designed to 
detect low-energy proton recoils, when, with 
large enough statistics, spectra of the type shown 
in Fig. 3 should reveal the presence of a pole 
operating in any accessible region of w*/?. 
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1*- PROTON SCATTERING CROSS SECTIONS BY CHEW-LOW EXTRAPOLATION* 


G. A. Smith,tt H. Courant, E. Fowler, H. Kraybill, J. Sandweiss, and H. Taft 
Physics Department, Yale University, New Haven, Connecticut 


Brookhaven National Laboratory, Upton, New York 
(Received November 17, 1960) 


The n’-p interaction cross section has been de- 
termined up to the 1-2 Bev region by many work- 
ers.'»? A new method recently suggested by 
Chew and Low’ allows one to infer the total 2* -p 
cross section at various energies by measuring 
the process p+p~-p+n+n"* at a specific labora- 
tory energy. In particular, consider the diagram 
of Fig. 1(a). In this diagram one of the protons 
emits a 1* meson which then scatters on the 
other proton; the pion emitter recoils as a neu- 
tron. The basic physical principle in the Chew- 
Low technique is related to the location and resi- 





q,, (P* me p,”)? 


due of a pole in the transition amplitude of the 
process p-n+z*. This pole is located in the 
nonphysical region of negative kinetic energy of 
the recoiling spectator neutron. The value of the 
residue at this pole is shown to be directly re- 
lated to the pion-proton cross section at a given 
pion-proton Q value. The value of this applica- 
tion of the Chew- Low technique to measure 1*-p 
cross sections is not in the actual measurement 
of the cross sections but rather as a test of the 


technique itself. 
Chew and Low suggest the construction of the 


function 


8 0( pnn*) 





. 
F(u?, p*) : an(w [i w* - $ w*(M ,*+ Le”) + 4(M,? = uF |* 


where w is the total energy of the proton and pion 
in their c.m. system and is related to their Q 
value by w=Q+M,+u,, p is the neutron recoil 
momentum, M, is the proton mass, M, is the 
neutron mass, J, is the 1* mass, qj, is the la- 
boratory momentum of the incident proton, and 


p,’ = -(M,/M,)[u.,? -(M, - M,)?). 


For any fixed value of w?, 870( pna*)/ap?aw? has 








re. 
7 n 
(a) 
FIG. 1. Feynman 
diagrams considered a 
in the Chew-Low anal- a" 
ysis of p+p—ptn+n*. -- 
P n 
« 
o 
P P 
(b) 


ap*8w* ”’ 





a second order pole at p=,” and the value of the 
residue at this pole is just F(w?,p,?). Thusa 
plot of F(w?, p”) vs p? extrapolated to p? = - yu,” 
(assuming M, =M,) will yield the residue which 

is related to the cross section in the following 
manner: 


F(w’?, = .”) = -2f%0 +,(@), 


where f is the pion-nucleon coupling constant 

and f?=0.08. It should also be noted that the 
diagram of Fig. 1(b) is also possible in 1* pro- 
duction. This process competes with that of Fig. 
1(a) and will provide background to the process 
of interest; however, this background should not 
affect the value of the residue F(w*, - u,”), where 
w=Q+M,+ yu, is for n* -p only. 

The 20-inch liquid hydrogen bubble chamber of 
the Brookhaven Cloud Chamber group was ex- 
posed to a beam of 2.85-Bev protons from the 
Cosmotron. 790 examples of p+p~p+n+n* have 
been identified in the analysis to date.* Each 
event was considered to be a possible candidate 
for the process of Fig. 1(a). Since the laboratory 
system and the rest system of the incident pro- 
ton are completely equivalent in the proton- 
proton collision, the event was treated in that 
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system in which the neutron momentum was 
lower. Since one must determine accurately the 
value of F(w?, p?) near p?=0 to perform the ex- 
trapolation into the unphysical region of p”, it is 
advantageous to select low neutron momenta. 

Because of the extreme reduction of statistics 
by the 8?0(pnz*)/ap?8w? factor of F(w*, p”), the 
final analysis was considered to be lacking suffi- 
cient statistics to measure more than two cross 
sections at different values of Q over rather wide 
intervals in p? and w*. The two intervals se- 
lected were for 62 <Q <226 Mev where the 1*-p 
cross section peaks in the 3/2, 3/2 resonance 
at 154 Mev to about 200 mb, and 226<Q < 793 
Mev where the 2*-p cross section falls from 
about 100 mb to 20 mb. The Chew-Low technique 
should measure some average cross section over 
each of these intervals. The results of this anal- 
ysis are shown in Fig. 2. 

Although the statistics make the extrapolation 
uncertain, on the basis of the least-squares 
linear and quadratic extrapolations shown in 
Fig. 2 the following conclusions seem to be justi- 
fied: (a) One can clearly detect the difference 
between the cross section for the region 62 <Q 
< 226 Mev and that which is found with the data 
in the region 226<Q <793 Mev, and in addition 
these cross sections have approximately the 
correct relative magnitude. (b) The data are 
consistent with the fact that 8?0( pnt) /ap? dw? 
must stay positive definite at least through 











p?= (??) min* 0.005 (Bev/c)?. 

So far as we know, this is the first Chew- Low 
extrapolation attempt, outside the deuteron case, 
where the answer has been known; in addition, 
it is felt that if a resonance exists in the 7,7 
case which is as large as the 3/2, 3/2 n*-p re- 
sonance, such an extrapolation should detect it. 
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knowledge the cooperation of Dr. Ralph P. Shutt 
and members of the Brookhaven Bubble Chamber 
group. 
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NEUTRINO EMISSION FROM BLACK-BODY RADIATION AT HIGH STELLAR TEMPERATURES 


Hong-Yee Chiu 
The Institute for Advanced Study, Princeton, New Jersey 


and 


P. Morrison 
Department of Physics and Newman Laboratory of Nuclear Studies, Cornell University, Ithaca, New York 
(Received November 28, 1960) 


That neutrino emission processes may affect 
the properties of matter at high temperatures, 
and hence stellar evolution, was first pointed 
out by Gamow and Schoenberg’ and by Ponte- 
corvo.” We list the known important neutrino 
emission processes as follows: 

(i) The so-called “urca” process, in which 
neutrinos are emitted by nuclear transitions.'»*>* 
This process depends only upon the observed 
nucleon-neutrino coupling, while all those which 
follow take place only with a universal Fermi 
interaction. 

(i) Neutrino bremsstrahlung”’®: 


e +(Z,A)-e +v+0+(Z,A). 


(iti) The photoneutrino or neutrino Compton 
effect*: 


y+e —e +V+D. 


We add two new processes, which are the sub- 
ject matter of this note®: 
(iv) Pair-annihilation neutrinos: 


e@ +e° —y+D. 

(v) Photon-photon neutrinos: 

l. y+y-v+D, 
2.y¥+y¥ ry tt. 


The energy lost from finite regions due to 
processes (i), (ii), and (é#i) is not important for 
temperatures below 10'°°K. For most densities, 
and for temperatures around 10°°K, processes 
(iv) and (v) are considerably more rapid. 

These last two processes involve direct neu- 
trino emission from black-body radiation itself, 
mediated either by virtual electron-positron 
pairs or by the real pairs present in thermal 
equilibrium. In spite of the more complex graphs 
which describe processes (v), they reach im- 
portance comparable with the simple process 
(iv) since the available energy in the center-of- 
mass system is typically higher, at least in the 
region of RT/mc*<«1. The graphs for these 
processes are shown in Fig. 1; process (v, 1) 














Vv. y, V y, - y 7. 
x ¥ (-] 7 , Va ‘ f 3 
“a a ~\ F ‘\ Ps 
\ fe G G Ja 
i Ja a Ja Ja 
e_ C+ Y % nm & 


yty—vt+y yty > ytv+y 
(1) (2) 
Processes (v) 


e+e*—> +7 
Process (iv) 


FIG. 1. Feynman diagrams for the processes con- 
sidered here. 


requires renormalization before it can be inter- 
preted. An estimate of process (v, 2) confirms 
the idea that it contributes to the energy loss 
rate about like process (iv) even though it is 
proportional to an additional three powers of 
1/137. 

Only the simple process (iv), with its unique 
Feynman diagram, has been calculated exactly. 
The transition amplitude is proportional to 
(ve)(@v) as in the theory of Feynman and Gell- 
Mann,’ from which we take the Fermi coupling 
constant, G. Calculation of the cross section 
for the pair process in the center-of-mass sys- 
tem is elementary, and gives 


=1.5x107*5 . 2 
o., 1.5x10 (c/y, Ep 1) cm’, 
where E7 is the total c.m. energy in units of 
m,c’, and v,¢, is the relative velocity in the c.m. 
system. The number of electron pairs present 
at equilibrium is well known.® In the pure vacuum 


radiation field, where the initial proper mass 
density is zero, the positron density becomes 


n ,Y =(1/n?)(/myc)f (8), B=(m,c?/kT) (1) 
and 


f (8) -{ ” exp[-A(1 +x?)¥*) x*dx 


0 
1 Ls] 
= 3) exp[-8 coshé]cosh26d@ 
oO 


= K,(8)/B. (2) 
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K,(8) is the modified Hankel function of second 
order.* The presence of nondegenerate matter 
expresses the density of positrons present at 
equilibrium, but does not affect the rate of pro- 
duction of positrons. The energy loss, however, 
depends only upon the production rate (and the 
equal annihilation rate), and is therefore inde- 
pendent of matter density, until degeneracy be- 
gins to affect the production of positron-electron 
pairs. For the nondegenerate case, then, the 
electron density does not matter; in the degen- 
erate case, if alson_>n,, we have 
(d) v : 

n =n, exp(-E _./RT); 
Ep is the electron Fermi energy, including 
mc. 

Note that 


(1 + $x?)>(1 +x?) >x. 


We find the approximations to (2): 


f (8) -e?( _— [-3x°B | x*dx 


=(1V2)(2/a)°/2eF, g>>1, 


fol8)={~ expl-x8]x*dx =2/8%, B«1. (3) 
0 


The rate of loss of energy density is given by the 
integral 


aU 
Pee E 3 3 

at fe. 2 ww" rei” $f 4 Py 
where the subscripts refer to electron and posi- 
tron quantities. Making such rough estimates 
as 


; = = 2 
E =E  =m,c ’ 


the integral can be approximated. We obtain 
-dU /dt = 7.5 x10*°[ f(8)]? ergs/cc-sec, 
(nondegenerate) 
=2.5x10'*(Z/A) exp(-E ,,/RT)f (8), 
(degenerate) 


where Z andA are mean values of atomic num- 
ber and atomic mass. 

In stellar problems the energy loss rate is 
usually expressed not in terms of unit volume, 
but per gram of matter, for which the relations 
above must be modified to read 


€ =(1/p)(dU /dt) ergs/g-sec. 


We list in Table I a few values of € for various 
processes of energy loss in the relevant tempera- 
ture and density region. It is to be observed that 
processes of type (v), so far only roughly esti- 
mated, are not only independent of density but 
even of degeneracy of the material. Both proc- 
esses (iv) and (v) are independent of chemical 
composition, apart from the effect of electron 
density on degeneracy. 

The thermal energy content of matter at T 
=10°°K is 10’’ ergs/gram. Therefore the relax- 
ation time for cooling by neutrino loss due to 
the processes we consider lies between 10° sec 
and 10*° sec. In any case, by the time T has 
reached 2 x10°°K, the relaxation time is only of 
the order of 10* to 10° sec. 

Temperatures of 2x10°°K or more are re- 
quired by the current theories of element syn- 
thesis for elements beyond iron, and by the 
ideas of supernova collapse.'°"* It seems that 
neutrino-loss processes change the means by 
which such temperatures may be reached, but 


Table I. Rates of energy loss by neutrino emission. . Energy loss in ergs/gram-second, at a density of 10° 











g/cc. 
Process 
Photon-photon 

Temperature Photoneutrino Annihilation (v ,2) teool 
(°K) (iii) (iv) (approximate) (sec) 

5x 108 102 10°? 103-5 1043 

1x 10° 104-2 107-4 10° 10° 

2x 10° 10° 1011-6 10" 10° 

2.5x 10° 1910-4 10!?-6 10” 104 








“In the last column is entered the relaxation time in seconds for losing the full thermal energy content of the 


material by neutrino emission. 
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do not modify their effects. 

The matter and thermal radiation in a stellar 
interior will remain in thermal equilibrium in 
spite of neutrino emission, provided that the 
energy loss rate is low enough to permit read- 
justment of the local near-equilibrium conditions. 
The virial theorem relates the-thermal energy 
strictly to the gravitational energy while this 
remains true. Neutrino loss will therefore not 
result in reduced temperature but only in the con- 
traction of the star, the gravitational energy 
making up for the neutrinos radiated. The central 
temperature will rise at a rate limited by neu- 
trino emission, until the time for gravitational 
readjustment —roughly measured by the free- 
fall times of 1 to 100 seconds (depending on 
density) —-becomes longer than the character - 
istic time of neutrino energy loss. Until central 
temperatures of about 3 x10°°K are reached, the 
star which has begun to emit neutrinos will every- 
where remain in a state of near equilibrium, and 
show increasingly rapid gravitational contrac- 
tion. This final stage of stellar evolution is 
thus controlled by the consequences of the uni- 
versal Fermi interaction, together with gravi- 
tational effects; nuclear processes do not play 
the determining role until the “phase change” 
from iron to helium is reached at 5x10°-7x10° °K. 

In view of the high neutrino emission rate, the 
first mark of supernovae might be high neutrino 
flux while the surface of the star has yet shown 
no sign of this sudden internal collapse. At 
T =5x10°°K, the neutrino flux from such stars 
will amount to 10°° neutrinos of average energy 
1 Mev per second. At a distance of 100 light 
years, the flux is around 10°*/cm?-sec. This is 
detectable.’* Therefore the establishment of a 


neutrino monitor station in terrestial or spatial 
laboratories may help us predict forthcoming 
supernovae. 

It is a great pleasure to thank Professors R. F. 
Christy, F. J. Dyson, W. A. Fowler, T. Gold, 
Mr. M. Levine, and Dr. T. T. Wu for much dis- 
cussion, and one of us thanks Professor J. R. 
Oppenheimer for the hospitality of the Institute 
for Advanced Study. 





‘G. Gamow and M. Schoenberg, Phys. Rev. 59, 539 
(1941). 

2B. M. Pontecorvo, J. Exptl. Theoret. Phys. 
(U.S.S.R.) 36, 1615 (1959) [translation: Soviet Phys. — 
JETP 9(36), 1148 (1959)]. 

3H. Y. Chiu, Ann. Phys. (to be published). 

4H. Y. Chiu and R. Stabler (to be published) . 

°G. M. Gandel’man and V. S. Pianev, J. Exptl. 
Theoret. Phys. (U.S.S.R.) 37, 1072 (1959) [translation: 
Soviet Phys. -JETP 10(37), 764 (1960)]. 

SCompare G. Wataghin, Phys. Rev. 66, 149 (1944) 
(before universal Fermi interaction). 

™R. P. Feynman and M, Gell-Mann, Phys. Rev. 109, 
193 (1958). 

8L. D. Landau and E. M. Lifshitz, Statistical Phys- 
ics, translation by R. F. Peierls and R. E. Peierls 
(Pergamon Press Ltd, London, 1958), p. 325. 

®*Tabulation of K,(8) can be found, e.g., in G. N. 
Watson, Theory of Bessel Functions (Macmillan Com- 
pany, New York, 1948), p. 737. 

105, A. Colgate and M. H. Johnson, University of 
California Radiation Laboratory Report UCRL-6097 
(unpublished) . 

‘lE, M. Burbidge, G. R. Burbidge, W. A. Fowler, 
and F. Hoyle, Revs. Modern Phys. 29, 547 (1957). 

“w, A. Fowler and F. Hoyle (to be published). 

134. G. W. Cameron, Astrophys. J. 130, 917 (1959). 
“PF. Reines and C. L. Cowan, Jr., Phys. Rev. 92, 
830 (1953). The neutrino background due to the sun is 

estimated to be less than 10'*/em?-sec on earth. 








575 








VoLUME 5, NUMBER 12 


PHYSICAL REVIEW LETTERS 


DECEMBER 15, 1960 





TOTAL CROSS SECTIONS OF PROTONS WITH MOMENTUM BETWEEN 10 AND 28 Gev/c 


A. Ashmore,* G. Cocconi, A. N. Diddens, and A. M. Wetherell 
CERN, Geneva, Switzerland 
(Received November 21, 1960) 


A beam of protons elastically scattered by the 
internal aluminum target of the CERN proton 
synchrotron! (scattering angle ~25 mrad) has 
been used for counter measurements of proton 
cross sections on protons, neutrons, and com- 
plex nuclei. 

The momentum of the proton beam was varied 
from 10 to 28 Gev/c by changing the operating 
energy of the accelerator and by isolating the 
elastic peak by means of magnetic selection. A 
typical momentum spectrum of the beam at 24- 
Gev operating energy is given in Fig. 1. It was 
obtained at ~80 m from the target, with a lm 
magnet and four scintillators defining the in- 
coming and outgoing beams (resolution of the 
system, 4p/p=5%). The spectra at other en- 
ergies were similar to this one. No particles 
were found around the elastic peak when the field 
of the analyzing magnet was reversed, thus con- 
firming that the particles detected at normal 
setting were indeed protons. 

The cross-section measurements were made 
utilizing the protons around (+ 2.5%) the elastic 
peak. At 80 m from the target their intensity 
was, at all energies, ~510° cm~? pulse", 
when 2 x10"! protons per pulse were circulating 
in the accelerator (+ pulse sec™’ at the highest 
energy, 1 sec™* at the lowest; duration of each 
pulse =~ 30 msec). 

The total cross sections on protons and neu- 
trons were measured with a “good geometry” 
setup. After magnetic selection, the incoming 
beam was defined by two scintillators, each 1 cm 
in diameter, 10 m apart, and the beam attenua- 
tion was measured by a row of 5 scintillators, 

2 cm in diameter, aligned with the first two and 
subtending, at the absorber, angles ranging 
from 2 to 10 mrad. The choice of these angles 
was dictated by the divergence of the beam, the 
expected multiple Coulomb scattering of the pro- 
tons in the absorber (~1 radiation length thick), 
and the expected distribution of the elastically 
scattered protons (the characteristic angle of 
the nucleon shadow scattering at 25 Gev is ~10 
mrad). 

For the p-p cross section, attenuation measure- 
ments were carried out alternatively with a car- 
bon and a polyethylene absorber. For the p-n 
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measurements H,O and D,O were used. Conven- 
tional coincidence circuits with resolving times 
of ~2x10-* sec measured simultaneously real 
and chance coincidences. 

The total cross section was obtained, at each 
energy, by extrapolating to zero solid angle the 
attenuations measured with the various aper- 
tures. For angles small in comparison with those 
expected for shadow scattering, this extrapola- 
tion should be linear in a solid angle plot, and in 
fact at all energies the experimental points after 
correction for multiple scattering? and beam di- 
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FIG. 1. Momentum spectrum of the external beam 


of protons elastically scattered in the target of the 
accelerator. Operating energy 24 Gev; Al target; 
scattering angle 25 mrad. The dashed line shows the 
shape expected for a monochromatic line of 24.2 
Gev/c. 
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FIG. 2. p-p and p-n 
total cross sections (dif- 
fraction scattering in- 
cluded) as a function of 
proton momentum. (Von 
Dardel et al., reference 
4, Longo et al., refer- 
ence 6.) 
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vergence did not deviate too much from a straight 
line. The slopes of these lines were, as ex- 
pected, inversely proportional to the square of 
the proton momentum.* 

The final results for the cross sections are 
given in Fig. 2, together with those of other au- 
thors. The correction for multiple Coulomb scat- 
tering and beam divergence was <44%, that for 
chance coincidence <4%. The errors quoted for 
our points are ~3 times larger than the purely 
statistical ones because of the contribution of the 
systematic errors introduced by these correc- 
tions and the extrapolation to zero solid angle. 
The relative values at different momenta should 
thus have errors smaller than those given in 
Fig. 2. Our values for the p-p cross sections are 
in good agreement, around 10 Gev/c, with those 
obtained by Von Dardel et al.* using a liquid hy- 
drogen absorber, and at 24 Gev/c with that ob- 
tained with the CERN 30-cm hydrogen bubble 
chamber.® The striking feature of the data con- 
tained in Fig. 2 is the constancy, at around 40 
mb, of the p-p cross section from 5 to 28 Gev/c. 

Taking into account the shadowing of the neu- 
tron by the accompanying proton in the deuteron, 
it seems that the p-n cross section is essentially 
the same as the p-p cross section at the three 
momenta where measurements were performed. 

At 24 Gev a “bad geometry” measurement of 
the cross section of protons on nuclei was per- 
formed. In this case only three scintillators were 


used after the absorber and subtended angles of 
15, 20, and 30 mrad. At these angles the effect 
of the shadow scattering was small and the other 
corrections were negligible. A linear extrapola- 
tion to zero solid angle based on these three 
measurements was taken as an estimate of the 
absorption (total minus diffraction scattering) 
cross section. The results of this measurement 
are displayed in Fig. 3. It is remarkable how 
well an A** law fits the data, though for the 
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FIG. 3. Nuclear absorption cross sections and 
mean free paths of 24.2-Gev/c protons in various 
elements. 
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lightest elements a good deal of transparency 
could be expected. 

We are extremely grateful to the Machine 
Group of the proton synchrotron for their very 
efficient operation of the accelerator at all en- 
ergies, and at all times. We thank Mr. C. A. 
Stahlbrandt for his assistance with the electron- 
ics. One of us (A.A.) wishes to thank CERN for 
making his visit possible. 
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An RESONANCE AND POSSIBLE DETERMINATION OF THE KYN PARITIES* 


Marc Rosst 
Universita di Roma, Rome, Italy 


and 


Gordon Shaw 
University of California, La Jolla, California 
(Received November 14, 1960) 


Experiments by Alston et al.’ and Ferro-Luzzi 
et al. ? at Berkeley, observing the reaction 


K~+p-A+nm7 +7", (1) 


have recently been reported which indicate a 
resonance in the A+z system at about 1380-Mev 
total energy, i.e., 50 Mev below the threshold 
for K+N. We propose that this resonance be 
described in terms of a quasi-bound state of the 
isotopic spin J=1, s-wave KN system. (The angu- 
lar distribution indicates J=4 for the resonant 
system, in agreement with this proposal.) Using 
our effective-range formalism® we will first 
establish the compatibility of this description 
with present data. Then we suggest how better 
data will permit fuller use of the effective-range 
theory to make detailed tests of this association 
of the Am resonance with the KN s state. Deter- 
mination of the orbital angular momentum of the 
resonant Az state (or the resonant 27 state that 
will also be present), and thus determination of 
the KAN (or KEN) parity, is discussed. 

This description can be decided by reference 
to the observed K” +p reactions. Let us consider 
process (1) at a fixed high energy and examine 
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those events where the 7*, say, has sufficiently 
high momentum so that it is definitely not asso- 
ciated with the A in this resonance. Such condi- 
tions were obtained in the experiment of Alston 
et al. We will speak in terms of the energy of the 
An” = system. We assume that the main production 
process near 1380 Mev involves two-body pro- 
duction of the 7* and the resonant state, with the 
subsequent decay of the latter. Call the ampli- 
tude for this entire process T,,,. In this picture, 
the quasi-bound KN system decays into A+7 

(or +7), the rapid energy dependence in T 
being associated with this decay. This depen- 
dence is characterized essentially by the two- 
body amplitude* 


[ky @7(K +N = A+n)| 


a >5 “a 
ht ’ (2) 


N 





1 y2 
=| im 3 /(1+ ais) 
a, 





with a,(ky) the J=1 KN complex scattering length, 
ky the KN momentum, extended below threshold 
where k—ix (x>0) and |a|*=|7T,(NZ)}?/|T(NA)|’. 
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Thus the rapid energy dependence is given by 
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The form (3) is explicitly that given if the 7” is 
considered to be “knocked off” the proton at first 
and to remain a “spectator.”* Equation (3) yields, 
by the way, an infinite derivative in the cross 
section in the form of a step (not a cusp) at the 
KN threshold (for Rea positive). This seems to 
agree with observation." 

Now we need to know the energy dependence of 
1/a,(x). Consider a system of coupled two-body 
channels in a given state of J, J,, I, parity, etc. 
Then the scattering amplitude can be written in 
the matrix form* 


al 1 
21+1 11+ 4h 


1 
T =k! * 2(y -ik (4) 


where M is a real symmetric (m xn) matrix which 
is an even function of all the channel momenta k; 
and k! +2 has entries (k/ + 2)5;= 6; h,!i* 2 (where 
below the ith threshold k;-ik;, «>0). It was 
shown in reference 3 that analogously to the one- 
channel effective-range theory the diagonal ele- 
ments are accurately given by an expression 
quadratic in k;: 


1 “aig, 2 
i 


. 2 
M ,{E) =M , {E 9) + 2¢ -k, (Eo) (5) 


1; 
with c,=1, c,=-3, etc., where R; is the range of 
forces in channel i and E, is any convenient en- 
ergy at which we make our expansion. In addition 
the nondiagonal elements M;; are to a good ap- 
proximation energy independent: M;,(E)=M ij(Eo)- 
The s-wave KN, J=1 state is coupled to two 
other channels (=z and Az) so that in order to 
use the full formalism we need to know 6 “zero- 
range” parameters M;;(E,). However, the rough 
low-energy data are only good enough to determine 
three J/=1 “zero-range” parameters: a, and the 
(Z/A);=1 ratio |a|*. Now as a first approxima- 
tion we take ks and k, to be constant over the 
energy region of interest. Then 


1 1 ’ 

= a oe! (6) 
a, (ky) a, NN 

where Ry is the range of forces in the KN chan- 


nel which we take to be LEN =0.4x107 cm. 
Among the scattering-length solutions of Dalitz 


and Tuan® only (a-) with 
a, =(1.0 +0.18%) x1078 em (7) 


could yield a (narrow) resonance for /=1. It is 
very important to note that only by looking at the 
energy dependence of a,(ky), e.g., noting from (6) 
that a,>Ry, can we decide that the positive a, 

[of the (a-) solution as contrasted to the (b-) 
solution] is due to an attraction strong enough to 
bind in the KN channel and not a repulsion or a 
strong nondiagonal interaction (which will yield 

a positive Rea but not a narrow resonance).*»* 
Inserting (6) and (7) into (3), we find that |T 
is peaked at 


|? 
res 


Ky =260 Mev/c, 


whereas experimentally ky = 180 Mev/c, and the 
full width at half maximum is 


Ax =120 Mev/c, 


whereas experimentally Ax=130 Mev/c. Some 
caution is called for in assessing this agreement 
since, besides our neglect of the energy depen- 
dence of ks and ka, (a) there is a direct or three- 
body amplitude interfering with T, es which may 
shift the peak position, although it should not dis- 
turb the width very much; (b) the scattering length 
(a-) used was based on too big an |a(ky =0)|?. We 
note that our proposal predicts only an J=1 reso- 
nance: Association of the resonance with the KN 
channel rules out an J=2 resonance and the scat- 
tering-length solutions (based on present data) 
permit resonances in either J=1 or 0 but not both. 
An important experiment needed to test our 
proposal in greater detail is the Coulomb inter- 
ference study on elastic K~ +p scattering: Con- 
structive interference would eliminate the scat- 
tering-length solutions with positive real part 
and invalidate our proposal. The present data’ 
are close to giving an answer. More precise low- 
energy K +p reaction experiments performed at 
several energies, e.g., accurate determination 
of the 2*x~/Z"n* ratio as a function of energy, 
would permit fuller use of the effective-range 
formalism (4) and (5). [A detailed survey is now 
in progress using this formalism to see just how 
much the effective-range terms can contribute 
to the K +p cross sections as well as the effects 
of the KYN parities; i.e., we are trying to de- 
termine how much we can tie down the quantities 
in (4) and (5) with a given accuracy of the K” +p 
data.] We note that a large energy region above 
threshold as well as the entire resonance region 
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below ky =0 is included in the range 
2)p2_p2 - < 
R; Ik, R, (ky O)i <1 


throughout which effective-range theory might be 
expected to be good.* 

From associating the Am resonance with the KN 
s state, it follows that if we can measure the 
orbital angular momentum of the Az (or the reso- 
nant £7 state) the KAN (or K=N) parity will be 
measured. One possibility of determining /, is 
from the A polarization. Referring to (1), we 
measure @ from the perpendicular to the plane 
of production of the resonant Az and the spectator 
particle, and consider s,,. and p,. Am systems. 
If P is the degree of polarization of the resonant 
systems, the A polarization is 


Si,: P, 


Piz: P(cos*@ - sin*6). 


Of course P depends on the production mechanism, 


angle, and energy: One may not readily find a 
situation with large P. 

A different possibility is to look at the energy 
dependence of |a|* over the resonance. The 27 
threshold is very favorably located close to the 
resonance, so that with any detailed data, s- or 
p-wave resonances could be distinguished. The 
present data already provide some indication of 
the KIN parity. Experimentally,’ |a|? is found 
to be roughly 1 at the KN threshold, whereas 
averaged over the resonant region, *|a|?<1. 
Theoretically, we have from (4),® in first approx- 


imation, 


lal?ak 


2l5 +1 21, +1 
: je, , (8) 


Thus, from (8), we find that |a|? decreases from 
its value at the KN threshold to that at the peak of 
the resonance by a factor of 2 for Jy =1, =1 (even 
KYN parities). On the other hand, if we con- 
sider />=0, |a|* remains essentially constant 

in this energy region, which may be in disagree- 
ment with experiment. 





*This work was performed in part under the auspices 
of the U. S. Atomic Energy Commission. 

tOn leave from Indiana University , Bloomington, 
Indiana. 
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UNIFIED APPROACH TO HIGH- AND LOW-ENERGY STRONG INTERACTIONS 
ON THE BASIS OF THE MANDELSTAM REPRESENTATION* 


G. F. Chew and S. C. Frautschi 
Lawrence Radiation Laboratory and Department of Physics, University of California, Berkeley, California 
(Received November 28, 1960) 


We wish here to outline an approach to the 
theory of strong interactions, based on the Man- 
delstam representation, that treats high- and 
low-energy phenomena in a unified manner. Our 
approach, which will be described in detail else- 
where, extends somewhat the original program 
proposed by Mandelstam.' The most striking 
achievement thus far is the inference, from the 
observed constancy of high-energy total cross 
sections, that in the low-energy elastic region 
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the P-wave interaction should be strong; at the 
same time, bound states or large low-energy 
phase shifts for J > 2 seem inconsistent with the 
theory. Furthermore, the divergence difficulty 
encountered by Chew and Mandelstam? in hand- 
ling low-energy P resonances appears to be re- 
moved by the Pomeranchuk relation® between 
high-energy particle and antiparticle cross sec- 
tions. 

We use the pion-pion interaction to illustrate 
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FIG. 1. The Mandelstam diagram 
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total energy—in pion mass units. s=0 
The assignment of locations to p;, 

P», and p3; corresponds to the ampli- 

tude A(s,t,u) of reference 4. Cor- 

responding assignments for B(s,t, u) 

and C(s,t,u) are obtained by inter- 

change of channel labels. ‘ 


our approach, although it will be clear that the 
essential features may be generalized. Figure 1 
shows the Mandelstam diagram for 7 elastic 
scattering, where the variables s, ¢, and u have 
the usual meaning.* The wedge-shaped physical 
regions are labeled as such, and the shaded 
areas indicate where the double spectral func- 
tions fail to vanish. Our central assumption is 
that the double spectral functions in the heavily 
shaded areas dominate those parts of the physi- 
cal regions whose distance from a boundary is 
of the order of magnitude of the width of the 
heavily shaded strips, i.e., 16 pion mass units 
(squared). This gives us a theory that covers 
not only low energies but arbitrarily high ener- 
gies with low momentum transfer. 

The above assumption seems a priori reason- 
able on the usual geometrical grounds, but in 
addition it is fortified by two empirical circum- 
stances. First, total cross sections generally 
are largest in the low-energy elastic region. 
Second, diffraction peaks in elastic scattering 
are always observed in the forward direction at 
high energies. The latter means, for example, 
that near the right-hand boundary of the upper 
physical region in Fig. 1 (the s region) the quan- 
tity, ImA(s,?#), for s large and fixed, falls off 
sharply as ¢ decreases from zero (forward di- 
rection). The width of the diffraction peak in NN 
and aN scattering corresponds to At~20, and it 
will be surprising if the same is not true for a7 
scattering as well as for other strongly interact- 


Physical 


Physical 











ing particle combinations. Now, we have 


1 (~.,, p(t’, s) 
ima(s,1)~> f at' ,s : 
1 t’-t 





4 


so a t dependence of the kind required for diffrac - 
tion implies a concentration of the double spec- 
tral function within a distance from the boundary 
of the order At. Such a concentration in the 
double spectral region labeled p, in Fig. 1, when 
applied to the ¢ channel (the lower right-hand 
physical region), is in agreement with the first 
of the above-mentioned empirical circumstances, 
i.e., the largeness of total cross sections at low 
energy. We shall return to this point at the end 
of this Letter. 

Fortunately, the double spectral functions in 
the heavily shaded areas are determined by re- 
latively tractable elastic unitarity conditions, as 
first pointed out by Mandelstam’ and recently 
re-emphasized by Cutkosky.* Because of space 
limitations we do not write these conditions here, 
although they will play a central role in actual 
dynamical calculations. Suffice it to say that in 
the mm problem there are three distinct “strip” 
functions to be calculated, which we have indi- 
cated in Fig. 1 as p,, p,, andp,. Each of these 
is determined exactly in the strip regions by an 
integral over bilinear combinations of absorptive 
parts of elastic amplitudes. It is only in the in- 
terior regions that inelastic amplitudes explicitly 
play a role. 
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In his original paper’ Mandelstam stated the 
opinion that attempting to calculate the parts of 
the double spectral function on which we are now 
focusing attention is not worthwhile because in- 
elastic effects are implied, in contradiction with 
the basic approximation, which is elastic. We 
believe that there is no inconsistency; the elastic 
unitarity condition is to be employed only where 
it is correct, in order to calculate double spec- 
tral functions in the strip regions. Inelastic 
scattering is not calculated completely, but that 
part occurring at low momentum transfers (im- 
plied, for example, in the s physical region in 
Fig. 1 by the existence of p,) should be well ap- 
proximated. 

It is worthwhile digressing momentarily to re 
late these ideas to those recently expressed by 
Salzman and Salzman* and by Drell,’ as well as 
by Pomeranchuk.® If we focus attention on the s 
physical region of Fig. 1, the double spectral 
function p, corresponds to diagrams in which 
only two particles are present in intermediate 
states but any number may be exchanged. In 
other words, p, is calculated from the Cutkosky 
diagram® shown in Fig. 2(a) and represents 
purely elastic effects in the s channel. On the 
other hand, p, is calculated from diagram 2(b), 
in which any number of particles are allowed in 
intermediate states but only two are exchanged 
(it is elastic in the ¢ channel). Obviously, then, 
we are calculating here the diffraction scatter - 
ing associated with inelastic transitions in which 
a single pion is exchanged. This is just the 
mechanism of Salzman and Salzman, Drell, and 
Pomeranchuk. We believe that our approach is 
more systematic, since it raises no questions it 
cannot answer about cross sections in unphysical 


(a) (b) 
FIG. 2. The two Cutkosky diagrams needed to cal- 


culate the double spectral functions in the strip re- 
gions. 
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regions; however, we can only discuss total 
cross sections and elastic scattering, not the 
distribution of inelastic events. 

Returning to the main argument, we now con- 
sider a second postulate: that total cross sec- 
tions for strongly interacting particles asymptot- 
ically approach constants at very high energies. 
Even within the scheme of approximation pro- 
posed here we believe that our combined require- 
ments of unitarity and analyticity are inconsist- 
ent with total cross sections that increase inde- 
finitely with energy. In any event, such a situa- 
tion seems nonsensical and would, we believe, 
be impossible if unitarity were completely en- 
forced. On the other hand, there are almost 
certainly solutions, such as the S-dominant solu- 
tions discussed by Chew, Mandelstam, and 
Noyes,” for which the asymptotic behavior of the 
double spectral functions implies a total cross 
section vanishing strongly at infinity. This al- 
ternative also we reject as unreasonable, feeling 
that a characteristic of strong interactions is a 
capacity to “saturate” the unitarity condition at 
high energies for those states that overcome the 
centrifugal barrier. Such a saturation should 
lead to a constant total cross section. 

We add a final related postulate, that the first 
diffraction peak approaches a constant limiting 
shape at very high energies, and thus we arrive 
at the general assumption”: 


lim ImA(s,/f-sf(t). (1) 


s-—> © 


-20<t<0 


From Fig. 1, this is equivalent to 


, 1 [ p,(t’,s) .,, , 
ca Lf mits) dt'~sf(t), (2) 
-20< t<0 


which seems unnatural unless 


lim p,(t,s)=sp,“(t). (3) 
s--@ 


4<t< 20 


Although limit (3) cannot be said to follow strictly 
from limit (1), we adopt (3) as the basis for fur- 
ther argument. Similar reasoning can be em- 
ployed for p,(t,s) and p,(t,s) to define p,/(t) and 
p,~(t). Then we borrow Pomeranchuk’s argu- 
ment® about consistency of the fixed ¢ dispersion 
relations to conclude that 


ps (t)=p,"(t). (4) 
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This equality corresponds to the statement that 
the total (7*z~) and (x*z*) cross sections become 
asymptotically the same. 

It is important to remark that the behaviors (3) 
and (4) appear to be consistent in the strip re- 
gions with the integral unitarity conditions which 
are to be used to calculate the double spectral 
functions. Thus we believe that we shall be able 
to find solutions of our equations satisfying these 
conditions. How many free parameters are pre- 
sent in such solutions we shall not know until 
further study has been made of the consistency 
of our equations in the interior regions where 
they are not exact and where a cutoff probably 
will be needed. 

We close by pointing out that condition (3), 
when the substitution law is invoked, implies 
that amplitudes at low energy are asymptotically 
proportional to tas t~«. Thus, regardless of 
how many parameters are allowed, S-dominant 
solutions are ruled out and a strong interaction 
in the P wave (and no higher waves) is implied." 
[At the same time, condition (4) appears to elim- 
inate the inconsistency encountered by Chew and 
Mandelstam? when they attempted to ignore the 
double spectral functions in a consideration of 
P-dominant solutions.] It seems, therefore, that 
one has in these considerations the beginning of 
an understanding of why strong interactions are 
so uniform in their manifestations: One may 
conjecture that the defining characteristic of 
strong interactions is that they are “as strong as 
possible.” The consequence from such a state- 
ment of a constant behavior for cross sections at 
high energies has always been plausible. Now 
we see that at low energies it appears inconsist- 
ent with general principles to have interactions 
so strong as to produce resonances or bound 
states ford >2. The “strongest possible” inter- 
action should, however, produce large low- 
energy P phase shifts. There remains an open 
question, of course, as to what extent one can 
determine the precise interaction strength on 


the basis of the “strip approximation” proposed 
here. 
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NEW APPROACH TO PERTURBATION THEORY IN QUANTUM ELECTRODYNAMICS 


Ivo Bialynicki-Birula* 
University of Rochester, Rochester, New York 
(Received November 4, 1960) 


It has been pointed out many times that the in- 
finite results in the perturbation approach to 
quantum field theory may be due to the singular 
character of the solution at zero value of the 
coupling constant. In this note we want to prove 
that this does happen. Our procedure is to calcu- 
late exactly the S matrix as far as this singular 
dependence on the coupling constant is concerned 
and to apply perturbation expansion to the re- 
maining nonsingular dependence. All results 
turn out to be finite’ in every order of perturba- 
tion theory as long as the coupling constant is 
kept finite. The measurable results do not de- 
pend critically on the coupling constant (and are 
the same as those obtained with the use of the 
ordinary renormalized perturbation theory), but 
the renormalization constants do. 

We will discuss our procedure in the case of 
quantum electrodynamics. Our starting point is 
the ordinary Lagrangian, 


£(x) = Flin mJ BA ANY -endy-Ay, (1) 


’ 


and we will not make any changes in the theory. 
Next we divide? this Lagrangian into the free part 
and the interaction part in a different way from 
the usual one, namely, 


£9 («) = Diy - m)y 
“ A~*J(iy-d . m)8y -3 aes (2) 
£ je) = -egy-Ay - (Sm) py 
+A~*Y(éy-8 - m)*%, (3) 


where A is a constant, whose role will become 
clear later. The free equations for the electron 
field can be solved exactly and they lead to a 
propagator having the Pauli-Villars structure, 


S(p) =(m - y-p)~* [1 - A-2(m - y-p)?] 3 


= (m - yp)? -3(m+h-y-p)* -3(m-A-y-p)"*. (4) 
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An indefinite metric appears in the world of vir- 
tual particles only and all observable (final and 
initial) states are characterized by a positive- 
definite metric in every order of perturbation 
theory. Let us now assume that 


AvsE “a. 


2(N +1) 
max 0 


where N is the highest order of perturbation 
theory, and E,,,,, is the highest energy of the 
initial and final particles that we are interested 
in. The S matrix can now be expanded into a 
power series in é,, and all terms are now finite 
provided the dependence of the electron propa- 
gator on é, is treated exactly. The additional 
term in the interaction Lagrangian can always 
be neglected in perturbation theory if the ener- 
gies of initial and final particles are finite, which 
means that this procedure cannot give answers 
concerning the limiting high-energy behavior. 
When the. energies of the initial or final particles 
tend to infinity the additional term in the interac- 
tion Lagrangian cannot be treated by perturba- 
tion theory. Ali renormalization constants ap- 
pear to have a singular dependence on the coup- 
ling constant of logarithmic type and they cannot 
be expanded into a power series ine,. All meas- 
urable quantities are given by a power series in 
é which is at least asymptotic. 

The procedure outlined above can be applied to 
all nonrenormalizable theories as well as to re- 
normalizable ones. A more detailed account will 
be published elsewhere. 





*On leave of absence from Warsaw University, 
Warsaw, Poland. 

'We will not discuss the infrared divergencies here. 

*There are, of course, many ways of dividing the 
initial Lagrangian which give the same answer. 
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INVESTIGATION OF THE FIVE-POINT FUNCTION* 


L. F. Cook, Jr.,? and Jan Tarskit 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received November 18, 1960) 


In this Letter we wish to describe and discuss 
the results of an investigation concerning the 
singularities of Feynman amplitudes for produc- 
tion processes. 

Our investigation was motivated by three con- 
siderations. First, experiments’ are being con- 
ducted in which five-particle processes are anal- 
yzed for the purpose of studying the 2-7 interac- 
tion,? and, in particular, the conjectured 1-7 re- 
sonance.* It is intended to extrapolate the data 
in at least one of these experiments (see Fig. 1), 
and a knowledge of the singularities of the ampli- 
tude is necessary to interpret properly the re- 
sults of such an extrapolation. Second, such 
production amplitudes are of interest in them- 
selves, and, also, are needed for the purpose of 
studying the three-particle intermediate states 
in dispersion theory. Third, Cutkosky has pointed 
out that the Feynman amplitude containing five ex- 
ternal momenta joined by a single internal loop 
has simple poles, which correspond to all the in- 
ternal momenta on the mass shell.‘ It is thus of 
interest to know whether these poles can be use- 
ful in theoretical analyses, and whether they 
have direct experimental effects; e.g., such 
poles could be near the pole indicated in Fig. 1.° 

In our investigation we restricted ourselves to 
single-loop diagrams without crossed lines, and 
in addition, we considered in detail the corre- 
sponding amplitudes for only four specific pro- 
cesses and for only one configuration of momenta 
for each process. Our reasons for considering 


FIG. 1. The diagram 
associated with a single- 
particle poleintheam- = \__ O 
plitude (N2x|Nx). An 
extrapolation to this pole \ 
is expected to give in- \ 
formation about the 1-7 \ 
interaction. \ 


only such diagrams are as follows: On the one 
hand, it can be easily shown that our general 
conclusions are not affected by the analytic pro- 
perties of the amplitudes which correspond to 
crossed diagrams, or to diagrams with a simpler 
internal structure. We may also note that the 
study of crossed diagrams does not lead to any 
new complications, and the simpler diagrams 
can be readily analyzed on the basis of previous 
studies.*»” On the other hand, diagrams with a 
more complicated internal structure yield singu- 
larities which are further removed from the re- 
gion of interest. Moreover, it has been shown 
that in some cases the essential features of such 
more complicated diagrams are also present in 
the single-loop diagrams for a given process,® 
and a similar conclusion may well be valid for 
production processes. 

We considered specific processes with specific 
configurations of momenta because a more gen- 
eral analysis seems to require a considerable 
amount of additional study. We felt that such a 
more general analysis could be more effectively 
carried out after a preliminary investigation has 
been made. 

The amplitudes which we considered in detail 
correspond to the diagrams of Fig. 2. Diagrams 
(a) and (b) are relevant to two of the experiments 
used for a study of the 7-7 interaction. Diagram 
(c) is pertinent for a possible experimental in- 
vestigation of the five-point poles. Since there 
are no single-particle poles in the amplitude 
{d2n\d7) analogous to that indicated in Fig. 1, 
it is likely that any behavior indicative of a pole 
is due to the five-point poles. Finally, diagram 
(d) is relevant in the study of the deuteron form 
factor.® We shall refer to the processes and the 
amplitudes of these diagrams as processes and 
amplitudes (a) - (d), respectively. 

We now consider a general Feynman amplitude, 
M,, corresponding to diagram (e) of Fig. 2, with 
momenta and masses as indicated. Such an am- 
plitude is a function of ten kinematical invariants, 
b;;*, i<j, where the pj; are defined by 
(1) 


Pag Pi ina to *9)-1,5° 


(Spin dependence is irrelevant for our discussion. 
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Me 


(d) (e) 


FIG. 2. Single-loop diagrams for production pro- 
cesses. Diagrams (a) -(d) correspond to the ampli- 
tudes which were studied in detail, and diagram (e) 
shows the labeling of the lines. In diagram (b) we as- 
sume that two of the internal pions are in resonance 
and behave as a single particle of mass 3.5m,. (This 
diagram has been suggested by R. J. Eden, Proceed- 
ings of the Tenth Annual Rochester Conference on 
High-Energy Nuclear Physics, 1960 (to be published)]. 


We also assume that each m; is different from 
zero.) It is also convenient to introduce quanti- 


ties xij defined by 


CG. 2 2 ; ; 
bi; =m. +m, +2m mx, for i<j, (2a) 
=X..5 x= 71. (2b) 
We form the matrix (x;;), and we are able to de- 
scribe concisely the possible locations of all the 
singularities of M, by the following theorem”: 
Every singularity of M, lies on a hypersurface 
(in the space of the complex invariants) defined 
by setting the determinant of a principal minor of 
the matrix (x;;) equal to zero. 

The various types of singularities, correspond- 
ing to the various hypersurfaces, can be givena 
simple physical interpretation: The order of a 
given principal minor equals the number of in- 
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ternal momenta on the mass shell. In particular, 
the five-point poles lie on the hypersurface 

det (x ;;) =0. In case of other singularities, there 
are internal momenta which do not affect the 
singularities. The singularities can then be 
identified with those of appropriate contracted 
diagrams. (A diagram is contracted by removing 
one or more internal lines and joining the end- 
points of each removed line. The result is a 
diagram with fewer vertices.) It should be em- 
phasized that every singularity lies on some such 
hypersurface, but the converse statement does 
not hold. In fact, the difficulty of the problem 
lies in determining which of these points are re- 
gular and which are singular on a suitably de- 
fined physical sheet. 

The complete details of our investigation will 
be submitted for publication in the near future, 
but a brief summary can be given here. For 
convenience we consider two of the invariants as 
(complex) variables, and the others as param- 
eters. We first choose such values for the param- 
eters that there is a region in the real plane of 
the two variables where the amplitude is real and 
analytic. For this choice of parameters the sin- 
gularities (for the two variables) can be readily 
obtained. We next vary the parameters to their 
physical values and follow the displacement of 
these singularities as the variation is done.®,’ 

We stated in the foregoing that we considered 
in detail only four processes and only one con- 
figuration of momenta for each process. We used 
the labeling shown in diagram (e) of Fig. 2 for the 
other diagrams, and we selected the following 
configuration in the center-of-mass system: 


Pi. = (E,2, q; 0, 0); - Psa = (Esa, 0, -3P, p), 
Pos = (Ex5, “qd; 0, 0); - Pas = (Eq; 0, -3P, -p), 
Py5= (E,5, 0, Pp, 0). (3) 


This configuration corresponds to momenta p,, 
and p,, incoming, and the others outgoing. Fur- 
thermore, we used the following incident labo- 
ratory energies: for process (a), T,=670 Mev; 
for process (b), 7, =648 Mev; for process (c), 
T,, = 650 Mev; and for process (d), T>= 500 Mev. 
[For processes (a) and (b), these energies are 
used in the experiments.] These incident ener- 
gies and Eqs. (3) now determine all the invariants. 
Some of our results are given in Fig. 3. This 
figure shows singularities in the complex p,,” 
plane, for each of the amplitudes (a) and (d). 
[Amplitudes (b) and (c) are described briefly in 
the sequel.] We chose the p,,” plane to display 
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FIG. 3. Singularities in the p,” planes associated, 
respectively, with diagrams (a) and (d) of Fig. 2. The 
values refer to the nucleon mass as a unit. See text 
for explanation of details. 


our results in order to compare the location of 
the five-point poles and other singularities of 
amplitude (a) with the pole indicate in Fig. 1, 
which occurs at p,,”=m,,”. For amplitude (d), 

the quantities of interest are the distances from 
the singularities to the region over which the am- 
plitude is integrated, since these distances are 
important in determining the effect of the singu- 
larities. The distances in the p,,” plane between 
the singularities and the location given by the 
value of p,,” in the physical process are illustra- 
tive, for a first orientation. The values of the 
other invariants are the actual physical values. 
However, we made the following exception: Pro- 
cess (d) is relevant for a study of the deuteron 
form factor, for which the amplitude (37|dd) is 
needed in an unphysical region, and with regard 
to this we somewhat arbitrarily reduced the 
value of p,,” to 10m,”. 

On the diagrams of Fig. 3 the crosses show the 
values of p,,” associated with the physical pro- 
cesses. The circle in diagram (a) indicates the 
pole implied by Fig. 1. Other geometrical fig- 
ures indicate in an obvious way the various types 
of singularities related to the single-loop dia- 
grams. In particular, the pentagons indicate the 
five-point poles, and the figures formed by two 
arcs indicate the thresholds for the production of 


two real nucleons; to these singularities corre- 
spond diagrams with lines 1, 3, and 4 contracted. 
Dotted lines indicate that a given pole lies on an 
unphysical sheet, and that a continuation along 
the indicated path will lead to the pole. (Other 
singularities on unphysical sheets are not dis- 
played.) Figures which are slightly removed 
from the real axis indicate that the points are 
singular in the limit Imp,,?~0, but the limit 
must be taken from the indicated half-plane. In 
this connection we may note that all the singula- 
rities other than the five-point poles (and the 
single-particle pole) are branch points, logarith- 
mic or square-root, and cuts must be prescribed 
which connect them and the point at infinity. Such 
branch cuts are shown on the diagrams of Fig. 3. 

The properties of singularities of amplitudes 
(b) and (c) do not differ greatly from those of the 
singularities of amplitude (a), and for this rea- 
son we do not give separate diagrams for these 
amplitudes. In particular, the five-point poles 
are located on an unphysical sheet, and in fact 
near the physical region, as in Fig. 3(a). How- 
ever, amplitudes (b) and (c) differ from ampli- 
tude (a) in having vertex-type singularities (anom- 
alous thresholds). Such singularities are present 
in amplitude (d), but not in amplitude (a). [We 
also note that the complete amplitude for process 
(b) has a single-particle pole at p,,?=my,’*. We 
recall that in the amplitude for process (c) there 
is no such pole. | 

The possible effects of a pole depend upon its 
residue as well as its distance from the relevant 
region, and we remark here that Cutkosky* has 
given a method of determining the residues of 
the five-point poles. A further discussion of the 
residues, however, is planned for a separate 
communication. 

We now turn to a brief discussion of our re- 
sults. Of immediate interest are the results for 
processes (a) and (b). Our investigation suggests 
very strongly that the five-point poles and the 
other singularities of amplitudes (a) and (b), 
respectively, have no noticeable direct effect on 
these processes. This conclusion is based on 
calculations for only special configurations of 
momenta, but even if moderate changes in these 
configurations are made, the essential aspects 
of the arrangement of the singularities should 
remain as in Fig. 3(a). Further, even if for 
some extreme configurations of momenta the 
amplitude should be markedly affected by the 
singularities, the contribution to the cross sec- 
tion of such configurations would not be great. 
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We also remark that the five-point poles have 
very little direct effect for process (c), for the 
chosen configuration of momenta. In contrast to 
the remarks in the preceding paragraph, it is of 
interest to determine whether certain extreme 
configurations exist for which the amplitude is 
noticeably affected by the five-point poles. Such 
effects could be experimentally measured, and 
the absence of a single-particle pole would sim- 
plify the theoretical interpretation. 

Finally, for process (d) (with the unphysical 
value p,,”= 10m”), we notice that there is a five- 
point pole on the physical sheet, and it can be 
expected to stay there for moderate changes of 
the configurations of momenta. Since the dis- 
tances from the singularities, including the pole, 
to the region of interest are such that their ef- 
fect on the amplitude can be significant, a fur- 
ther study of this process is suggested. On the 
other hand, a preliminary calculation suggests 
that the amplitude (37 |NN) has no singularities 
other than the normal thresholds near the region 
of interest; in particular, the relevant five-point 
poles lie on an unphysical sheet. Thus it appears 
that these other (anomalous) singularities can be 
ignored in the study of the three-pion contribu- 
tion to the nucleon form factor. 

In conclusion, we should like to make one gen- 
eral observation which may be of interest. 
Studies of the vertex and scattering amplitudes® 
have shown that the presence of external particles 
which are near instability leads to singularities 
other than the normal thresholds. This pattern 
of behavior is also illustrated by our results re- 
garding the amplitudes relevant in the nucleon 
and deuteron form factors, viz., in the latter 
case the anomalous singularities are on the phys- 
ical sheet. 
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at the Lawrence Radiation Laboratory, where 
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fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 








PERTURBATION THEORY FOR AN INFINITE 
MEDIUM OF FERMIONS. I. Abraham Klein, 
University of Pennsylvania, Philadelphia, Penn- 
sylvania (Received October 13, 1960). 


New forms of the linked-cluster expansion for 
the ground-state energy of a normal spherical 
system are given in terms of the exact single- 
particle propagator and of the associated self- 
energy operator. These forms are shown to be 
stationary with respect to variations of the latter 
quantities. A new concise proof of the Hugenholtz- 
Van Hove theorem on single-particle energies is 
then possible. The variational principle also ex- 
plains the insensitivity of recent binding energy 
calculations to the choice of single-particle en- 
ergies. It is shown how all results may be ex- 
tended to systems with anisotropic features. 


PERTURBATION THEORY FOR AN INFINITE 
MEDIUM OF FERMIONS. II. DERIVATION OF 
THE LANDAU THEORY OF FERMI LIQUIDS. 
Abraham Klein, University of Pennsylvania, 
Philadelphia, Pennsylvania (Received October 
13, 1960). 


The methods employed in the previous papers 
of this series for the ground-state energy of a 
normal Fermion system are extended to finite 
temperatures. The basic assumption is that the 
density matrix of the actual system in thermo- 
dynamic equilibrium can be derived by adiabatic 
transformation from the density matrix of a sys- 
tem of noninteracting particles at the same tem- 
perature, volume, and chemical potential. How- 
ever, the spectrum of the noninteracting system 
must be dependent on the intensive parameters. 
The self-consistency of the scheme is exhibited 
by demonstrating that the effective potential char - 
acterizing this system can be related to the self- 
energy operator of the actual system. Two equi- 
valent versions of the theory are derived. In one 


of these the basic equations are precisely those 
of the Landau theory of Fermi liquids with the 
difference, however, that all quantities are given 
explicitly from first principles. 


COULOMB SCATTERING OF POLARIZED ELEC- 
TRONS. W. R. Johnson, T. A. Weber, andC. J. 
Mullin, University of Notre Dame, Notre Dame, 
Indiana (Received August 19, 1960). 


An approximate analytical expression for the 
spin-dependent cross section for Coulomb scat- 
tering, valid to order (aZ)*, is developed. The 
development is based on a modification of the 
Sommerfeld-Maue wave function. The (a Z)? and 
(a Z)*® terms in the resulting cross section are 
identical with the corresponding terms found by 
the second Born approximation. The (a Z)* term 
is verified by analytically summing the corres- 
ponding term in the Mott series. Graphs, com- 
paring the approximate analytical expressions 
for both the cross section and asymmetry function 
with exact numerical results, are included. 


CYCLOTRON RADIATION FROM RELATIVISTIC 
PARTICLES WITH AN ARBITRARY VELOCITY 
DISTRIBUTION. Ludwig Oster, Laboratory of 
Marine Physics, Yale University, New Haven, 
Connecticut (Received October 5, 1960). 


A formula is derived which represents the 
spectral and angular behavior of cyclotron ra- 
diation emitted by a relativistic particle moving 
in a constant magnetic field and having arbitrary 
velocity components parallel and perpendicular 
to the field. This formula is used to describe 
the emission of an assembly of particles having 
a Maxwellian distribution function. All broaden- 
ing mechanisms of practical importance-—colli- 
sions, relativistic mass variability, and Doppler 
effects—are included and discussed. 


ANALYTIC PROPERTIES OF SINGLE-PARTICLE 
PROPAGATORS FOR MANY-FERMION SYSTEMS. 
J. M. Luttinger, Columbia University, New York, 
New York (Received October 11, 1960). 


Certain general properties of single-particle 
propagators for a system of interacting fermions 
are derived. In addition, the properties of the 
proper self-energy part G,(¢) which were used 
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in previous work on the ground-state energy and 
on the Fermi surface are established. In parti- 
cular, the fact that to all orders of perturbation 
theory in the interaction, ImG,(x -i0*) behaves 
like Cp(x - )? (Cp>0) for x very near yp, is 
proved. 


EFFECT OF ELECTRON EXCHANGE ON THE 
DISPERSION RELATION OF PLASMONS. Oldwig 
von Roos and Jonas S. Zmuidzinas, Jet Propul- 
sion Laboratory, California Institute of Techno- 
logy, Pasadena, California (Received October 5, 
1960). 


A calculation of the influence of electron ex- 
change on the dispersion relation of a high-density 
electron gas at 0°K is described. The result is 
compared with those obtained by various authors 
using different methods. 


CORRECTION TO THE DEBYE-HUCKEL 
THEORY. John G. Trulio and Stephen G. Brush, 
Lawrence Radiation Laboratory, University of 
California, Livermore, California (Received 
October 13, 1960). 


The correction to the equation of state of a 
classical electron gas is calculated by the method 
of Abe. The results are compared with those 
recently obtained by a different method by Bowers 
and Salpeter. 


SOLUTION OF THE FUNCTIONAL DIFFEREN- 
TIAL EQUATION FOR THE STATISTICAL EQUI- 
LIBRIUM OF A CRYSTAL. Robert M. Lewis and 
Joseph B. Keller, Institute of Mathematical Sci- 
ences, New York University, New York, New 
York (Received October 13, 1960). 


The s-particle distribution functions (s = 1, 2, +++) 
of classical equilibrium statistical mechanics are 
determined for a crystal, as power series in the 
temperature. They are obtained by solving Bogo- 
lyubov’s functional differential equation. From 
the distribution functions, the thermodynamic 
functions of a crystal are computed as power 
series in the temperature. The leading terms 
in these series are the usual classical results 
which are customarily derived by assuming that 
the potential energy is a quadratic function of the 
particle displacements. The further terms, which 
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depend upon the nonquadratic or anharmonic terms 
in the potential, provide corrections to the usual 
results, which become more important as the 
temperature increases. If only a few terms in 
the series are used, the results will be valid at 
temperatures low compared to some characteris- 
tic temperature of the crystal, e.g., the melting 
temperature. Since they are based on classical 
mechanics, the results are valid only at tempera- 
tures high compared to the Debye temperature. 

The series expansions of the distribution func- 
tions and thermodynamic functions may be viewed 
as the low-temperature analogs of the virial ex- 
pansions, which are low-density expansions. As 
in the case of the virial expansions, all the terms 
are determined explicitly in analytic form, but 
their actual evaluation is difficult. 


DEFECTS IN IRRADIATED SILICON. I. ELEC- 
TRON SPIN RESONANCE OF THE Si-A CENTER. 
G. D. Watkins and J. W. Corbett, General Elec- 
tric Research Laboratory, Schenectady, New 
York (Received October 7, 1960). 


The Si-A center is a major, radiation-damage 
defect produced in “pulled” silicon by a room tem- 
perature irradiation. As a result of studies de- 
scribed in this paper (I), and the following one 
(II), it is concluded that this center is a lattice 
vacancy with an oxygen atom impurity bridging 
two of the four broken bonds associated with the 
vacancy. Spin resonance and electrical activity 
arise from an electron trapped in the other two 
bonds. In this paper (I), the spin resonance 
studies are described. A molecular orbital treat- 
ment of the trapped electron wave function satis- 
factorily accounts for the observed g tensor, as 
well as the hyperfine interaction observed with 
neighboring 4.7% abundant Si®® nuclei. Study of 
the changes in the spectrum of a sample subjected 
to uniaxial stress are also described. Under 
stress, the amplitudes of the individual resonance 
components (which correspond to different orien- 
tations of the defect in the crystal) are observed 
to change. This results from (1) electronic re- 
distribution of the trapped electrons among the 
defects, and (2) thermally activated reorientation 
of the defects themselves under the applied stress. 
These two effects are separated and a quantitative 
study of their magnitudes and signs, as well as 
their rates, is given. The results confirm many 
of the important microscopic features of the model. 
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DEFECTS IN IRRADIATED SILICON. II. IN- 
FRARED ABSORPTION OF THE Si-A CENTER. 
J. W. Corbett, G. D. Watkins, R. M. Chrenko, 
and R. S. McDonald, General Electric Research 
Laboratory, Schenectady, New York (Received 
October 7, 1960). 


The Si-A center is a major, radiation-damage 
defect produced in “pulled” silicon by room tem- 
perature irradiation. In this paper (II), we pre- 
sent the infrared measurements which, in con- 
junction with the spin resonance measurements 
of the preceding paper (I), establish the identity 
of the Si-A center. A new infrared absorption 
band is observed at 12 yu in electron-irradiated 
Silicon. This band is shown to be a vibrational 
band of impurity oxygen in the lattice. Macro- 
scopic and microscopic correlations between the 
12u band and the spin resonance of the Si-A 
center are presented. The macroscopic corre- 
lations are of production rate, recovery, etc. 
The microscopic correlations derive from the 
absorption of polarized infrared radiation by sam- 
ples of various crystallographic orientations, sub- 
jected to a uniaxial, compressive stress. Partial 
alignment of the defects is induced by the stress 
and is detected as a dichroism in the 12y band. 
This alignment is compared to the corresponding 
alignment studies in the spin-resonance meas- 
urements in paper I. It is shown that the kinetics 
and magnitude of the response to the stress are 
the same for the defects observed in both types 
of measurements. This shows that the 12y band 
arises from the Si-A center and establishes the 
configuration of the oxygen in the defect. These 
results, together with the results of paper I, 
allow us to conclude that the Si-A center is a 
lattice vacancy with an oxygen atom bridging two 
of the four broken bonds associated with the va- 
cancy. The remaining two bonds can trap an elec- 
tron, giving rise to the spin resonance spectrum 
of the defect. The identification of the Si-A 
center indicates that the vacancy is mobile ina 
room-temperature irradiation. 


ELECTRICAL PROPERTIES OF n-TYPE GaSb. 
Alan J. Strauss, Lincoln Laboratory, Massa- 
chusetts Institute of Technology, Lexington, 
Massachusetts (Received October 3, 1960). 


The electrical conductivity and Hall coefficient 
of n-type GaSb doped with Se or Te have been 
measured, in most cases at 77°K and 300°K, for 
samples with net donor concentrations between 


about 6 x10" cm™® and 2x10"* cm~’. In general, 
the data are consistent with the two-band model 
which Sagar has proposed for the conduction 
band of GaSb, but systematic differences are 
observed between the properties of Se-doped 
and Te-doped samples. It seems likely that _ 
these differences are associated with impurity 
conduction of the metallic type. For the Te- 
doped samples, the electron mobility at 77°K 
apparently varies in an anomalous manner with 
increasing impurity concentration. 


EFFECT OF STATIC STRAINS ON DIFFUSION. 
L. A. Girifalco and H. H. Grimes, Lewis Re- 
search Center, National Aeronautics and Space 
Administration, Cleveland, Ohio (Received 
September 16, 1960; revised manuscript re- 
ceived November 7, 1960). 


A theory is developed that gives the diffusion 
coefficient in strained systems as an exponential 
function of the strain. This theory starts with 
the statistical theory of the atomic jump frequency 
as developed by Vineyard. The parameter deter - 
mining the effect of strain on diffusion is related 
to the changes in the interatomic forces with 
strain. Comparison of the theory with published 
experimental results for the effect of pressure 
on diffusion shows that the experiments agree 
with the form of the theoretical equation in all 
cases within experimental error. 


ULTRASONIC ATTENUATION IN NORMAL 
METALS AT LOW TEMPERATURES. A. B. 
Bhatia and R. A. Moore, University of Alberta, 
Edmonton, Canada (Received September 6, 
1960; revised manuscript received October 24, 
1960). 


Expressions for the attenuation a, and a; of 
plane dilatational and shear sound waves are 
obtained by solving the Boltzmann transport 
equation for the electron distribution function f 
without assuming the existence of a relaxation 
time 7 for the collision term in this equation. 
Instead the collision integral is considered to 
arise explicitly from the interaction of electrons 
with thermal phonons and impurities. Making 
the usual “ideal metal” assumptions, it is found 
that the attenuation in general depends on a set 
of effective relaxation times 7 LM Which are 
associated with the various terms in the expan- 
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sion of f in a series of spherical harmonics 

Y tm; ); the ty yy are independent of the sub- 
script M, and hence the same set {77} deter- 
mines both ag and a;. Explicit expressions for 
T, are derived. 

‘por the case of all 7 L equal to one another 
and equal to 7, say, the analytical expressions 
for a, and a; obtained here are the same as 
those of Pippard. However, usually the 7; are 
not equal to one another. It is then found that 
when >> / (A is the wavelength of the sound 
wave and / a mean free path of the electrons), 

7 in Pippard’s expressions must be replaced by 
T, and, contrary to what is usually assumed, a 
would not be in general proportional to the elec- 
trical conductivity 0 (oa 7,). When r<<l, the 
attenuation, with one exception, is independent 
of {7 L and is the same as that given by Pippard. 
For \~l, and the Tz not equal to one another, 

a@ may be calculated numerically if the ratios 
T,/T, are known; the results of one such calcu- 
lation show that the deviations from Pippard’s 
analytical expressions are at most about 20%, 
provided 7 in the latter is identified as 1,. 

Lastly, the possible influence of electron- 
electron collisions on attenuation is briefly dis- 
cussed. 


DOMAIN CONFIGURATIONS ABOUT NONMAG- 
NETIC PARTICLES IN IRON. William D. Nix 
and Robert A. Huggins, Department of Materials 
Science, Stanford University, Stanford, Cali- 
fornia (Received October 3, 1960). 


The closure domain configuration about a non- 
magnetic particle in iron is considered in terms 
of magnetostatic, interfacial wall, and magneto- 


strictive energies. Equations for the total energy 


of the domain configuration and the value of the 
geometrical parameter @ have been derived. The 
effects of particle size and shape are considered. 


BLEACHING AND RECOVERY OF F CENTERS 
IN KCl. Loretta E. Silverman* and L. I. Gross- 
weiner, Department of Physics, Illinois Institute 
of Technology, Technology Center, Chicago, 
Illinois (Received October 12, 1960). 


F centers in KCl are optically bleached at room 


temperature with a minimum light energy of 20.1 
+1.0 ev per F center. The F band is recovered 
on further room-temperature x-raying with a 
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minimum energy of 69+ 5 ev per F center. The 
increase in the M and R bands during optical 
bleaching of the F band is nullified upon the re- 
covery of the F band by further x-raying. The 
recovery kinetics are second-order in the number 
of optically bleached F centers with a rate-con- 
stant of (7.94 0.8)x10-2° cm’/F center-sec. The 
rate law requires an energy transfer from the 
bulk lattice and is consistent with an exciton me- 
chanism. 


*Present address: Electronics Research Division, 
Armour Research Foundation, Technology Center, 
Chicago, Illinois. 


CROSS-RELAXATION AND CONCENTRATION 
EFFECTS IN RUBY. Roy W. Roberts, Melabs, 
Palo Alto, California, James H. Burgess, Stan- 
ford University, Palo Alto, California, and 
Harold D. Tenney, Melabs, Palo Alto, California 
(Received October 13, 1960). 


Cross relaxation effects in ruby maser crystals 
are treated by introduction of a cross- relaxation 
probability in the rate equations. Detailed solu- 
tions have been obtained for several specific pro- 
cesses and compared to recent experiments. It 
is shown that cross relaxation can improve maser 
performance even in the absence of impurity dop- 
ing. Pulse experiments at 0.06% and 0.14% chro- 
mium ion concentrations in a ruby traveling-wave 
maser are interpreted in terms of a five-spin 
process in addition to a four-spin process. 


PHOTOCONDUCTIVITY OF SILVER CHLORIDE 
CRYSTALS UNDER PULSED X-RAY IRRADIA- 
TION. Alwin E. Michel,* Department of Engi- 
neering Physics, Cornell University, Ithaca, 
New York (Received October 5, 1960). 


The photoconductivity produced by irradiation 
with 0.2-sec x-ray pulses was studied in AgCl 
as a function of temperature (80 - 280°K), x-ray 
intensity and penetration, field strength, and 
crystal preparation, in order to obtain informa- 
tion about lifetimes and mobilities of electrons 
and holes. The measurements of electron life- 
times (~1yusec) and mobilities in air-grown crys- 
tals are in agreement with those reported in the 
literature. The photoconductive response can be 
described by assuming deep electron traps through- 
out the volume of the crystal and possibly a dis- 
turbed surface layer. No hole motion is observed 
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below 250°K; above that temperature the schubweg 
per unit field is estimated at 5x10-® cm?/v. 

The electron lifetime in crystals grown and 
annealed in He is much smaller than in the air 
grown samples. Assuming the same mobility in 
both samples, the lifetime at 80°K is 3x10~° sec. 
At higher temperatures the pulses show long tails, 
and between 200 and 280°K the saturation time 
varies exponentially with 1/7. It is assumed that 
there exist shallow traps (~ 0.08 ev) in a thin sur- 
face region which otherwise has a long electron 
lifetime as compared to the bulk of the crystal. 
Measurements on crystals doped with 20 ppm of 
Cut indicate that the Cut ions do not act as effec- 
tive electron traps. On the other hand, the pres- 
ence of 1 ppm of Ni ions reduces the lifetime at 
30°K to less than 3x10°"! sec, indicating a cap- 


ture cross section of the Ni ion larger than 300 A?. 


*Present address: IBM Research Laboratory, Pough- 
keepsie, New York. ‘ 


DEGENERACY IN Ag,Te. C. Wood, V. Harrap, 
and W. M. Kane, Research Division, Philco 
Corporation, Philadelphia, Pennsylvania (Re- 
ceived August 24, 1960; revised manuscript 
received November 8, 1960). 


The Hall coefficient, resistivity, and Seebeck 
coefficient of n- and p-type specimens of Ag,Te 
have been measured over the temperature range 
from 55 to 300°K. These results indicate that 
the compound is highly degenerate over the whole 
temperature range studied. Calculations were 
made of the effective masses, mobility ratios, 
and energy gap and gave order-of- magnitude 
values. 


SURFACE STRUCTURES AND PROPERTIES OF 
DIAMOND-STRUCTURE SEMICONDUCTORS. 

D. Haneman,* Barus Physics Laboratory, Brown 
University, Providence, Rhode Island (Received 
October 6, 1960). 


Low-energy electron-diffraction and secondary 
electron emission measurements have been made 
on (111) and (111) surfaces of GaSb and (100) sur- 
faces of InSb. To account for the diffraction pat- 
terns observed both for these materials and pre- 
viously for Ge and Si, a general model for (111) 
surfaces of diamond-structure semiconductors 
is proposed. Every second atom, counting along 
alternate close-spaced rows, is raised with re- 


spect to its neighbors, being bonded to the sub- 
layer by three p bonds while the “dangling bond” 
is s-type. The remaining three-fourths of the 
surface atoms have dangling p bonds and are 
bonded to the sublayer by trigonal sp?-type bonds. 

The Ga or (111) face of GaSb has maximum | 
sticking coefficients of 10-° and 10~* for oxygen 
and CO,, respectively, these values being 10 
times greater than those found for the Sb or (111) 
face. Multilayer adsorption of oxygen takes place 
on all the surfaces measured. The oxygen can 
be removed by heat treatment alone. Evidence 
is presented to show that diffusion of oxygen into 
the bulk is an important mechanism for regenerat- 
ing the clean surfaces by heat treatment. Carbon 
dioxide adsorbs on GaSb so as to show structure. 
It apparently deposits as an unbroken molecule, 
and, unlike oxygen, does not build up several 
layers. 

*Now at Physics Department, University of New 
South Wales, Sydney, Australia. 


FRIEDEL SUM RULE FOR A SYSTEM OF IN- 
TERACTING ELECTRONS. J. S. Langer, 
Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania, and V. Ambegaokar,* Westing- 
house Research Laboratories, Pittsburgh, 
Pennsylvania (Received October 3, 1960). 


The Friedel sum rule is derived for a system 
of interacting electrons in a periodic potential. 


*Present address: Institute for Theoretical Physics, 
Copenhagen, Denmark. 


ABSORPTION SPECTRA OF F, , Cl, , Br,”, AND 
I, IN THE ALKALI HALIDES. C. J. Delbecq, 

W. Hayes,* and P. H. Yuster, Argonne National 
Laboratory, Argonne, Illinois (Received October 
12, 1960). 


After x-ray irradiation at liquid nitrogen tem- 
perature, holes are trapped by forming F,” mole- 
cule ions in LiF; Br,” molecule ions in KBr con- 
taining one of the impurities Ag*, Tl*, Pbt*, or 
NO,°; and L, molecule ions in KI containing one 
of the impurities listed above. An investigation 
of the optical and paramagnetic resonance spectra 
of crystals containing F,°, Br,°, and I, centers 
leads to an identification of the optical transitions 
of these centers. The polarizations and relative 
intensities of the optical absorptions in the X,~ 
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molecule ion series are discussed with reference 
to the energy level scheme to be expected for 
such a species. 


*Now at the Clarendon Laboratory, Oxford, England. 


SURFACE MOBILITY OF COPPER IONS ON 
CUPROUS OXIDE. Rudolf Frerichs and Irving 
Liberman, Northwestern University, Evanston, 
Illinois (Received September 9, 1960; revised 
manuscript received November 7, 1960). 


This paper reports on the motion of Cut ion va- 
cancies on the surface of cuprous oxide at room 
temperature under the application of an electric 
field. The measurement of the mobility of the Cut 
vacancies was made by means of a “time of flight” 
procedure. The formation of luminescent centers 
is the unique property of Cut vacancies that makes 
them directly observable. The mobility of the Cut 
vacancies at room temperature is about 107"! 
cm?/volt-sec. The variation of the mobility with 
temperatures between 28°C and 55°C is observed. 
From these data the constants of the diffusion 
equation D =D, exp(4H/RT) are computed: D, 
=5x10~-? cm? sec, AH = 8100 calories. The low val- 
ues obtained for these constants show that the 
ionic current follows low-resistance paths formed 
by the crystal grain boundaries or along the sur- 
face of the crystal. 


TUNNELING CURRENT IN ESAKI DIODES. 
Clayton Wilson Bates, Jr., Solid-State Physics 
Group, AVCO Research and Advanced Develop- 
ment Division, Wilmington, Massachusetts (Re- 
ceived October 4, 1960). 


The integral giving the net tunneling current 
flowing across the junction in an Esaki diode, 
E 


v 
I=A f {f,E)-f,(E)} Zp, (2)p, (ak, 
E 
c 


is evaluated under the normal assumptions that 
(é. -E,,) and (E,, - £) are of the order of 2kT. 
The resulting expression is 


qV/kT) 


qv /kT) . 


ANE, -E Yl -¢ 
a/2 





I=- 


(m+n )e +(l+e 


where A” is an arbitrary constant and m,n, and 
a are functions of the Fermi levels on both 
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sides of the junction, the location of the band 
edges, and the absolute temperature. This ex- 
pression is plotted as a function of the applied 
voltage for temperatures of 200°K, 300°K, and 
350°K for donor and acceptor concentrations of 
10° cm=* and 1.6x10'® cm~’, respectively. The 
resulting curves compare quite favorably with 
those of Esaki. 


ELECTRIC FIELD GRADIENTS IN DILUTE AL- 
LOYS OF SILVER. C. A. Giffels,* G. W. Hinman, 
and S. H. Vosko, f Department of Physics, Car- 
negie Institute of Technology, Pittsburgh, Penn- 
sylvania (Received August 31, 1960). 


The directional correlation of the gamma rays 
emitted in the decay of Cd*"* has been measured 
with the cadmium embedded in the cubic silver 
lattice in order to observe the electric field 
gradients produced by various solute atoms. The 
alloys studied contained small concentrations of 
cadmium, indium, tin, antimony, or germanium. 
The results were compared to various theories 
of screening and were found to be in agreement 
with a cos(2kr+)/r* falloff of shielding charge, 
where & is the wave number at the Fermi surface. 
This work corroborates recent results of Rowland 
by an independent method. 


*Now at Bell Telephone Laboratories, Whippany, 
New Jersey. 

+tNow at McMaster University, Hamilton, Ontario, 
Canada. 


EXCITONS AND PLASMONS IN SUPERCONDUC - 
TORS. A. Bardasis and J. R. Schrieffer, Uni- 
versity of Illinois, Urbana, Illinois (Received 
October 13, 1960). 


The Anderson-Rickayzen equations of motion 
for a superconductor derived within the random 
phase approximation (RPA) are used to investi - 
gate the collective excitations of superconductors. 
A spherical harmonic expansion is made of the 
two-body interaction potential V(k, kK’), and a 
spectrum of excitations whose energies lie within 
the energy gap 24 is obtained. These excitations 
may be characterized by the quantum numbers 
L and M involved in the potential expansion. For 
an L-state exciton to exist, the L-wave part of 
the potential must be attractive at the Fermi sur- 
face. Odd-L excitons have unit spin and may be 
considered as spin waves. For s-state pairing 
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in the superconducting ground state, the plasmon 
mode corresponds to the L = 0 exciton whose en- 
ergy is strongly modified by the long-range Cou- 
lomb interaction. For a general potential several 
bound states may exist for given L and M. If the 
L-wave potential is stronger than the s-wave 
part of the potential, the system is unstable with 
respect to formation of L-state excitons. In this 
case, the ground state is formed with L-state 
pairing, special cases of which are the p-state 
pairing considered by Fisher and the d-state 
pairing proposed recently by several authors 

for the ground state of He* and nuclear matter. 
Corrections to the Anderson-Rickayzen equations 
are discussed which lead to a new set of exciton 
states if the L-wave potential is repulsive. These 
excitons are interpreted as bound electron-hole 
pairs, as opposed to the particle-particle exci- 
tons present with an attractive L-wave potential. 


OPTICAL CONSTANTS OF SILVER. E. A. Taft 
and H. R. Philipp, General Electric Research 
Laboratory, Schenectady, New York (Received 
October 10, 1960). 


The optical constants for electrolytically pol- 
ished silver samples were determined from 
normal-incidence reflectance data in the spectral 
region 2 to 10 ev. The sharp minimum in the 
curve for k near 3.8 ev, which is associated with 
the onset of interband transitions, is deeper than 
in previous results. A second minimum is ob- 
served near 9.2 ev, which is approximately the 
calculated energy for simple free-electron plasma 
oscillations. 


OPTICAL-MODEL ANALYSIS OF LOW-ENERGY 
ELECTRON-HYDROGEN SCATTERING. B. A. 
Lippmann and H. M. Schey, Lawrence Radiation 
Laboratory, University of California, Livermore, 
California (Received September 30, 1960). 


The electron-hydrogen scattering problem is 
reduced to an equivalent one-body problem by 
the use of an optical-model potential. Two effects 
are examined quantitatively: the influence of the 
Pauli principle and the elastic scattering in the 
triplet state. For a long-range optical potential, 
the effect of the Pauli principle is small; however, 
for a short-range optical potential, the Pauli 
principle changes the zero-energy cross section 
by an order of magnitude. 


The elastic scattering results, though in agree- 
ment with many previous calculations, are dis- 
appointing. It is argued that this is due to an in- 
consistent application of the optical model, and 
that a more consistent use of the optical model 
is equivalent to (and therefore provides a deriva- 
tion for) the heuristic procedure, used by Martin, 
Seaton, and Wallace, which yields scattering 
lengths close to the Rosenberg, Spruch, and 
O’Malley bounds. 


IONIZATION CROSS SECTIONS FOR PROTONS 
ON HYDROGEN GAS IN THE ENERGY RANGE 
0.15 TO 1.10 Mev. J. W. Hooper, Engineering 
Experiment Station and School of Electrical En- 
gineering, Georgia Institute of Technology, 
Atlanta, Georgia and E. W. McDaniel, D. W. 
Martin, and D. S. Harmer, Engineering Experi- 
ment Station and School of Physics, Georgia In- 
stitute of Technology, Atlanta, Georgia (Received 
September 23, 1960). 


Measurements have been made of the ioniza- 
tion cross section for protons incident on hydro- 
gen gas in the energy range 0.15- 1.10 Mev. The 
experimental cross section in this region can be 
represented by 0; = (3.454 0.20)E~ (0.870 + 0.005) 

x 10°!” cm?/molecule, where E is the incident 
proton energy in Mev. The experimental results 
are in excellent agreement with a Born approxi- 
mation calculation, which is discussed. 


DETECTION OF DOUBLE RESONANCE BY 
FREQUENCY CHANGE: APPLICATION TO Hg™ 
Robert H. Kohler, Department of Physics and 
Research Laboratory of Electronics, Massachu- 
setts Institute of Technology, Cambridge, Massa- 
chusetts (Received September 27, 1960). 


A new type of double-resonance experiment that 
depends on wavelength effects rather than on 
polarization effects is discussed. Incident polar - 
ized light is replaced by incident light of wave- 
length coincident with just one component of the 
structure to be measured. The analyzer is re- 
placed by a cell of gas that absorbs just that 
same component and lets the others pass. Mag- 
netic resonance from the excited component to 
one of the others is monitored by increases in 
the light transmitted through the absorbing gas. 
This experiment requires that the Doppler width 
be smaller than the structure. 
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This method was first applied to measure the 
hyperfine structure of the *P, state of Hg*”. The 
incident light and absorbing gas were both sup- 
plied by separated Hg’®, whose resonance line 
coincides naturally with one component of the 
Hg” hyperfine structure. Measurement of Hg*™ 
is discussed in detail. These hfs intervals were 
found: f($— 3) =7551.613 + 0.013 Mc/sec, and 
F(¥— 3) =13 986.557 + 0.008 Mc/sec. The magnetic 
dipole and electric dipole interaction constants, 
calculated without quadratic hfs corrections, are 
a = -5454.569+ 0.003 Mc/sec, and b= -280.107 
+ 0.005 Mc/sec. 

Means for applying the method when there is 
no isotope coincidence are given. This new 
method is compared with the polarization tech- 
nique and is found to give signal-to-noise ratios 
that are orders of magnitude higher. 


MICROWAVE SPECTRUM OF CHLORINE DI- 
OXIDE. I. ROTATIONAL ASSIGNMENT. R. F. 
Curl, Jr., and James L. Kinsey,* Department 

of Chemistry, Rice University, Houston, Texas, 
John G. Baker,’ Department of Physical Chemis- 
try, Cambridge University, Cambridge, England, 
James C. Baird, + George R. Bird, and Robert F. 
Heidelberg, Department of Chemistry, Rice Uni- 
versity, Houston, Texas, and T. M. Sugden, D. R. 
Jenkins, and C. N. Kenney, Department of Physi- 
cal Chemistry, Cambridge University, Cambridge, 
England (Received October 13, 1960). 


The microwave spectrum of chlorine dioxide has 
been investigated and several rotational transitions 
assigned. The resulting rotational constants (in 
Mc/sec) are: 


C10,: A-=52072, 
CI"0,: A=50725, 


These constants are limited in precision by the 
presence of hyperfine structure which has been 
only approximately treated. The structural para- 
meters of chlorine dioxide as obtained from the 
rotational constants are 77)_9= 1.473+0.01 A; 
Z0-Cl-O=117° 36’+1°. Approximate hyperfine 
constants are given. 


C =8332; 
C = 8295. 


B=9952, 
B=9952, 


*Present address: Department of Chemistry, Uni- 
versity of California, Berkeley, California. 


+Present address: Bristol University, Bristol, England. 


t Present address: California Research Corporation, 
Richmond, California. 
|| Present address: Polaroid Corporation, Cambridge, 


Massachusetts. 
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a* CORRECTIONS TO HYPERFINE STRUCTURE 
IN HYDROGENIC ATOMS. Daniel E. Zwanziger,* 
Columbia University, New York, New York (Re- 

ceived August 26, 1960). 


The a* term in the ratio of the hyperfine split- 
ting in the 2S state of the one-electron atom to 
the hyperfine splitting in the 1S state is recalcu- 
lated, and a new theoretical value for this ratio 
is obtained which is in agreement with the ex- 
perimental value, thereby eliminating a pre- 
viously reported discrepancy. The calculation 
consists in the evaluation of the low-momentum 
parts, of order a*hfs, of the expression for the 
lowest order radiative level shift in the bound 
interaction representation with external Coulomb 
and magnetic dipole fields. By rearranging the 
terms so as to display the gauge invariance of 
the matrix elements with respect to the external 
potentials, considerable simplicity is achieved, 
and the formulas are easily interpreted as a 
generalization of the expression for the lowest 
order Lamb shift. The contribution from soft 
photon intermediate states is obtained by an ex- 
tension of the method developed by Schwartz 
and Tiemann for evaluating the Bethe logarithm, 
and an Appendix contains a tabulation of twelve 
analogous integrals which were integrated nu- 
merically, and which may be of use elsewhere. 
The calculated value of the ratio is 41.000 034 5 
+ 0.000 000 2) which agrees with the experimental 
values for hydrogen: 4(1.000034 6+ 0.000 000 3), 
and deuterium: #(1.000 034 2+ 0.000 000 6). 


*Now at the University of California, Berkeley, 
California. 


LOW-ENERGY PROTONS PRODUCED IN THE 
DEUTERON BOMBARDMENT OF NUCLEI. 

E. W. Hamburger,* B. L. Cohen, and R. E. Price, 
University of Pittsburgh, Pittsburgh, Pennsyl- 
vania (Received July 5, 1960). 


When nuclei are bombarded by fast deuterons 
a continuous spectrum of protons is observed at 
energies corresponding to Q<-2.2 Mev. The 
shape and angular distribution of this continuum 
were studied for several targets, from Li to Au, 
with 14.8-Mev deuterons. The angular distribu- 
tions are strongly peaked at or close to 0° for 
light elements; for heavy elements the peak 
broadens and moves to ~50°. The energy spectrum 
shows a peak which, with increasing angle, moves 
to lower energies for light elements, and moves 
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JRE to higher energies for heavy elements. All fea- the well-known prompt comparison method. 
er,* tures of the results are explained using a semi- Systematic errors can be reduced and greater 
2e- classical theory in which the deuteron is broken accuracy obtained for nuclei which exhibit both 
up in the external field of the nucleus. The a prompt and a delayed event. A transistorized 
: breakup occurs at or near the nuclear surface time-to-amplitude converter was used. The 
lit- for light elements, but quite far out from it in mean life of the second excited state of Ca**[(4.8 
, heavy elements. Observed cross sections are + 0.3) x10-!° sec] and the first excited state of 
“ail far larger than theoretical predictions for either Sc*’(< 5x10°"* sec) were measured. Also the 
” Coulomb or nuclear breakup. mean life of the first excited state of Hg'® [(3.5 
*Now at Departamento de Fisica, Faculdade de Filo- + 0.5) x 10°" sec] was determined by the usual 
" sofia, Ciencias e Letras, Universidade de Sao Paulo, self-comparison method. 
Sao Paulo, Brazil. *Present address: Department of Physics, Rutgers 
“A University, New Brunswick, New Jersey. 
mb NEUTRON-NEUTRON SCATTERING LENGTH. 
e Kirk W. McVoy,* Lawrence Radiation Laboratory, LONGITUDINAL POLARIZATION OF O* POSI- 
University of California, Livermore, California, TRONS. J. C. Hopkins,* J. B. Gerhart, F. H. 
ial and Brandeis University, Waltham, Massachu- Schmidt, and J. E. Stroth,t Department of Physics, 
d, setts (Received August 29, 1960). University of Washington, Seattle, Washington 
(Received October 3, 1960), 
t The final-state interaction of the two neutrons 
from the reaction 1~ +d ~-2n+y has a pronounced The longitudinal polarization of positrons from 
™ and distinctive effect on the momentum spectrum the pure Fermi transition of O* has been studied 
of the outgoing particles. In particular, the by the method of Bhabha scattering. A prismatic 
= neutron spectrum is sharply peaked in the neigh- beta-ray spectrometer selected 1.0-Mev positrons 
" ; borhood of 2 Mev, with a shape that is quite sen- which then were focussed by a solenoidal magnetic 
sitive to the strength of the mm interaction. In lens on a Supermendur foil magnetized by the lens. 
. the region of this peak, the relative neutron- Scattered positrons and recoil electrons were de- 
5 neutron momentum is so small that the mn inter- tected in fast coincidence. The scattering asym- 
tal action is completely characterized by its scat- metry observed upon reversal of the lens magnetic 
) tering length. Hence it is proposed that a meas- field was (4.164 0.80) % which corresponds to posi- 
‘ urement of the shape of the neutron spectrum tron longitudinal polarization +(0.97+4 0.19)v/c. 
from this reaction may provide a convenient *Now at Los Alamos Scientific Laboratory, Los 
means of measuring the neutron-neutron scatter - Alamos, New Mexico. 
ing length. Neutron spectra are calculated in an tNow at San Jose State College, San Jose, California. 
impulse approximation, for several assumed 
values of the scattering length. “t appears from 
their shapes that, in this way, present neutron- ASYMMETRIC FISSION OF BISMUTH. T. T. 
detection techniques should be capable of de- Sugihara and J. Roesmer,* Jeppson Laboratory 
‘ice, termining the scattering length to within 25 %. of Chemistry, Clark University, Worcester, 
= *Present address: Brandeis University, Waltham, Massachusetts, and J. W. Meadows, Jr.,f Cyclo- 
Massachusetts. tron Laboratory, Harvard University, Cambridge, 
Massachusetts (Received August 24, 1960; re- 
4 vised manuscript received November 7, 1960). 
NEW APPLICATION OF DELAYED COINCI- 
| DENCE TECHNIQUES FOR MEASURING LIFE- fu, *eYmmesie modo of mage tviion fn 
TIMES OF EXCITED NUCLEAR STATES: Ca 200 
\u, a " sion of Bi*”” with 36-Mev protons. About 0.3% of 
u- AND Sc". P. C. Simms, HN. Bencser-Koller, the fissions contribute to this mode. At 58 Mev 
and C. S. Wu, Columbia University, New York, e 
New York (Received October 11, 1960). no evidence for asymmetric fission (< 0.05% of 
. total fissions) as a separate mode could be found. 
trum A new application of the delayed coincidence The fission cross sections at 36 and 58 Mev are 
oves technique for measuring short nuclear lifetimes 1.9 and 11.3 mb, respectively. The narrowness 
eS has been developed which is basically similar to of the 36-Mev asymmetric peak leads to the sug- 
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gestion that the asymmetric fission of bismuth re- 
sults from the fission of a single nuclear species 
and from a closed-shell effect, similar to the 
fine structure observed in low-energy fission of 
heavy elements. This asymmetric fission is con- 
sidered to occur from states of relatively high 
excitation energy. However, the possibility of 
asymmetric fission also occurring from states 
of low excitation energy, whether following neu- 
tron evaporation or as a consequence of an in- 
elastic proton interaction, cannot be ruled out. 
The symmetric fission observed with both 36- 
and 58- Mev protons is consistent with the results 
obtained by Fairhall in the fission of bismuth with 
22- Mev deuterons. 

*Present address: Nuclear Science and Engineering 
Corporation, Pittsburgh, Pennsylvania. 

+Present address: Argonne National Laboratory, 
Argonne, Illinois. 


NEUTRON TRANSFER TO EXCITED STATES 

IN N** IN THE REACTION N“*(N"*, N*5)N?5. 
Kenneth S. Toth, Oak Ridge National Laboratory, 
Oak Ridge, Tennessee (Received September 9, 
1960). 


The neutron-transfer reaction, N**(N™*, N*°)N?°, 
was investigated from 6° to 60° c.m. with 28-Mev 
N“* ions accelerated in the Oak Ridge 63-inch 
cyclotron. Circular strips of aluminum foil, 
each encompassing a known angular increment, 
were used to stop the radioactive N** (10- min). 
Excited states in N** are unstable with respect 
to proton emission; the first excited state in 
N'® is 5.28 Mev above ground. Therefore, N** 
nuclei resulting from transfers to the ground 
state of N’® are at least 5 Mev more energetic 
than those resulting from transfers to N** ex- 
cited states. The two groups of N* particles 
were distinguished by placing suitable absorbers 
in front of the aluminum catcher foils. It was 
found that transfers to excited states become 
more abundant relative to ground-state transfers 
as the incident N“* energy is lowered. 


BETA-GAMMA DIRECTIONAL CORRE LATIONS 
IN THE DECAY Sb’ —- Te'**. Helmut Paul,” 
Physics Department, Purdue University, Lafa- 
yette, Indiana (Received October 10, 1960). 


We have measured the £-y directional correla- 
tions involving the 1.6-Mev 8 group of Sb'™ and 
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the second excited 2* state (at 1325 kev) of Te?*. 
Both of the y rays (0.72 and 1.32 Mev) de-exciting 
this state were used. For the 0.72-Mev y ray, the 
results for the coefficient € in W() =1+€P,(cos@) 
are: €=0.20+0.02, 0.18+0.02, 0.19+0.03, and 

0. 22+ 0.06 for 8 energies of 1.02, 1.16, 1.30, 

and 1.44+ 0.07 Mev, respectively. For the 1.32- 
Mev y ray, we find € = -0.28+ 0.08 for a B-energy 
range of 1.0-1.6 Mev. From these results, the 
mixing ratio for the 0.72-Mev y ray is 5 =+0.879-3. 
It is shown that the 1.6-Mev £8 transition must be 
due, at least partly, to the B;; matrix element. 
Assuming u =x =0 (in Kotani’s notation), our data 
can be fitted by either Y /z =0.04 (i.e., an es- 
sentially unique transition), or Y /z =1.2. 


+ 
Present address: Studiengesellschaft fiir Atomen- 
ergie, Lenaugasse 10, Wien VIII, Austria. 


LIQUID-DROP NUCLEAR MODEL WITH HIGH 
ANGULAR MOMENTUM. Robert Beringer and 
W. J. Knox, Department of Physics, Yale Uni- 
versity, New Haven, Connecticut (Received 
October 10, 1960). 


The equilibrium shapes of rotating liquid drops 
are calculated in the spheroidal approximation. 
The solutions are prolate forms for high angular 
momenta such as are produced in compound nu- 
clei resulting from heavy-ion bombardments. 
This is thought to account for the angular and en- 
ergy distributions observed in heavy-ion-induced 
fission and light-particle evaporation. 


MEASUREMENT AND INTERPRETATION OF 
NEUTRON TOTAL CROSS SECTIONS OF CAR- 
BON, CALCIUM, AND LEAD. R. M. Wilenzick, 
G. E. Mitchell,* K. K. Seth, and H. W. Lewis, 
Duke University, Durham, North Carolina (Re- 
ceived August 29, 1960). 


The neutron total cross sections of carbon, 
calcium, and lead have been determined in the 
energy region 200-960 kev using a time-of-flight 
technique. The observed cross section for car- 
bon has been analyzed in terms of the effective- 
range theory as well as the bound-level contri- 
bution with results which agree very well with 
those obtained by the (d,p) stripping experi- 
ments. Also it is concluded that the upper iimit 
for the total width of the “proposed” resonance 
at 610 kev is 100 ev. For calcium and lead the 
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data are analyzed for resonance parameters. It 
is found that the s-wave strength function for 
Ca*° increases rapidly with increasing energy. 
It is concluded also that the course of the Pb 
cross sections can be explained only by the pre- 
sence of a broad (#,*= 0.05) s-wave resonance 
in Pb, at E,,=5154 15 kev with resonance 
width ©, = 100+ 15 kev. 


*Now at Florida State University, Tallahassee, 
Florida. 


MEASUREMENT OF DEUTERON POLARIZA- 
TION PRODUCED BY d-a SCATTERING AT 

1.07 Mev. Lee G. Pondrom* and J. W. Daughtry, 
Aeronautical Research Laboratories, Wright- 
Patterson Air Force Base, Ohio (Received Oc- 
tober 13, 1960; revised manuscript received 
November 10, 1960). 


Elastic scattering from He* gas in the energy 
region of the 1.07-Mev resonance has been used 
to polarize deuterons. After scattering through 
30° in the laboratory by the polarizer, the polar- 
ization of the deuteron beam has been analyzed 
by the reaction He*(d, p)He* near 400 kev. This 
was accomplished by observing protons at 0° and 
90°, and measuring the dependence of the count- 
ing rate ratio on the deuteron polarization. 
He*(d, p)He* is well described in this energy re- 
gion by a single-level Breit-Wigner formula, 
and therefore has a predictable sensitivity to 
deuteron polarization. The experimental results 
are consistent with the magnitude of the polariza- 
tion component (7,,) calculated from the d-a 
phase-shift analysis. 


*Now at Physics Department, Columbia University, 
New York, New York. 


ANALYSIS OF K-MESON PRODUCTION BY p 
ANNIHILATION. T. F. Hoang,* Lawrence Ra- 
diation Laboratory, University of California, 
Berkeley, California (Received August 22, 1960; 
revised manuscript received November 17, 
1960). 


The production of K and 7 mesons in p-p anni- 
hilations is compared for antiprotons of 1-Bev/c 
momentum and antiprotons at 480-Mev kinetic 
energy. By modifying the Fermi model by as- 
suming a 7-7 resonant state with a mass slightly 
greater than two pion masses, good agreement 
with the results is obtained and the interaction 


volume 2, is no longer too large. An isobaric 
state with a mass greater than or equal to three 
pion masses is ruled out: 


*Now at Laboratoire de Physique, Ecole Polytech- 
nique, Paris, France. 


ANTIPROTON ANNIHILATIONS IN PROPANE. 
Sulamith Goldhaber, Gerson Goldhaber, Wilson 
M. Powell, and Rein Silberberg, Lawrence 
Radiation Laboratory and Department of Physics, 
University of California, Berkeley, California 
(Received September 1, 1960). 


An experiment to study the p annihilation pro- 
cess at 1.05 Bev/c was performed with the Law- 
rence Radiation Laboratory 30-in. propane bubble 
chamber. It was observed that the K-meson pro- 
duction in annihilation events rises sharply with 
the increase in energy, namely from (4+1)% for 
annihilations at or near “rest” to (8+1)%. On the 
other hand, the pion multiplicity was not observed 
to increase appreciably with the increase of avail- 
able energy. We have found a pion multiplicity of 
5.0+0.2. These numbers are discussed in this 
paper and compared with existing models for the 
p annihilation process. It is pointed out that with 
further increase in bombarding energy, different 
models may differ appreciably in the above quan- 
tities. 

We have observed a p -H annihilation cross 
section of 51+ 10 mb and a pf -C annihilation cross 
section of 368+ 60 mb at 1.05 Bev/c p momentum. 
Crude determinations of the p charge-exchange 
process ~—which turns out to be forward peaked— 
and of pf inelastic-scattering events leading to 
pion production are also discussed. 


HEAVY NUCLEI IN THE PRIMARY COSMIC 
RADIATION AT PRINCE ALBERT, CANADA. 

Il. H. Aizu,* Y. Fujimoto, S. Hasegawa, M. 
Koshiba,t I. Mito,} J. Nishimura, and K. Yokoi, 
Institute for Nuclear Study, University of Tokyo, 
Tanashi, Tokyo, Japan, and Marcel Schein, |! De- 
partment of Physics, Enrico Fermi Institute for 
Nuclear Studies, University of Chicago, Chicago, 
Illinois (Received August 19, 1960; revised 
manuscript received November 4, 1960). 


The investigation of the low-energy primary 
cosmic radiation has been extended to include 
the heavier elements of Z 29; the light elements 
Li, Be, and B; anda particles. The results of 
a previous paper on carbon, nitrogen, and oxy- 
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gen have also been confirmed with better statis- 
tics. The energy spectra of all these components 
show a general similarity in shape. A possible 
deviation of the light-element spectrum from this 
similarity is discussed. The abundances of var- 
ious elements in the low-energy region of 200 to 
700 Mev per nucleon are essentially the same as 
observed in the higher energy region. 


*Department of Physics, Rikkyo University, Tokyo, 
Japan; now at Department of Physics, Northwestern 
University, Evanston, Ilinois. 

tNow at Department of Physics, University of Chi- 
cago, Chicago, Illinois, on leave of absence from 
University of Tokyo, Tokyo, Japan. 

tcollege of General Education, University of Tokyo, 
Komaba, Tokyo, Japan. 

Deceased. 


DEMONSTRATION OF QUANTUM MECHANICS 
IN THE LARGE. T. B. Day, University of Mary- 
land, College Park, Maryland (Received Sep- 
tember 22, 1960). 


An example is given which demonstrates in a 
straightforward and dramatic manner that when 
two particles like (K°+K°) or (2y) are created 
simultaneously, the probabilities involved in ob- 
serving any further events related to their simul- 
taneous creation must be calculated quantum- 
mechanically and are correlated, even for macro- 
scopic distances in absorbing media. In parti- 
cular, a correlation in the polarizations of the 
two y rays from positronium annihilation as a 
function of the thickness of magnetized iron 
through which they are passed is pointed out by 
way of a proposed experiment. 


SHELL-MODEL CALCULATIONS OF THE 
HYPERFINE EFFECT IN ».-MESON CAPTURE. 
H. Uberall,* Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania (Received October 14, 
1960). 


Numerical values for the difference in u -meson 
capture rates from the two hyperfine states of a 
mesonic atom are obtained for F’®, Al’, and P™, 
employing the simple Mayer-Jensen version of 
the shell model for Al*’ and P™, and a more 
accurate version for F’*. In all these cases, we 
get considerably larger values than those found 
in previous estimates, thus rendering experi- 
ments more feasible. In view of the consequences 
of atomic conversion on the observation of the 
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capture rate differences, conversion effects in 
the rate and asymmetry of u-meson decay elec- 
trons are discussed also. 


*Present address: Harrison M. Randall Laboratory 
of Physics, University of Michigan, Ann Arbor, Michi- 
gan. 


MODEL OF HYPERON DECAY. L. Wolfenstein, 
Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania (Received October 7, 1960). 


The vanishing of the asymmetries in the decays 
= -n+n~ and =*—n+n* are explained in a model 
in which all = decays occur via virtual K-meson 
decays. The model then predicts a A-decay 
asymmetry of similar magnitude but opposite 
sign to that of 5*—p+7°. Certain other predic- 
tions of this model are discussed. 


INELASTIC SCATTERING OF A > HYPERON 
WITH AN EMULSION NUCLEUS. H. E. Fisk and 
D. J. Prowse, Physics Department, University 
of California, Los Angeles, California (Received 
October 3, 1960). 


During a systematic study of fast hyperons re- 
sulting from the nuclear capture of K” mesons, 
an event has been found which is interpreted as 
the elastic scattering of a =~ hyperon by a bound 
neutron. The reaction kinematics are insufficient 
to determine the nature of the nuclear potential 
for =~ hyperons but it appears that when about 
10 of these events have been found and analyzed, 
it will be possible to infer the sign of the potential 
from the general effect it will have on the observ- 
able kinematics of the events. 


MOMENTUM SPECTRUM OF MUONS. W. Pak,* 
Ss. Ozaki, f B. P. Roe,! and K. Greisen, Labora- 
tory of Nuclear Studies, Cornell University, 
Ithaca, New York (Received October 12, 1960). 


The momentum spectrum has been measured 
with a magnetic spectrometer for cosmic-ray 
particles incident vertically and at 68° zenith 
angle, in the range 2 to 70 Bev/c. The apparatus 
discriminates strongly against all particles but 
muons. The vertical muon spectrum is found to 
be in good agreement with the results reported 
by Pine et al. and by Ashton et al. No compara- 
ble data at 68° have been published. Assuming 





1960 


in 
aC - 


ory 
ichi- 


“ays 
10del 
son 


and 


ived 


re- 
Sy 


und 
‘ient 
al 


ed, 
ntial 
erv- 





VOLUME 5, NUMBER 12 


PHYSICAL REVIEW LETTERS 


DECEMBER 15, 1960 





muon creation to occur entirely through pion 
decay, a single-pion production spectrum ac- 
counts for the muon spectra at both zenith angles. 
The muon positive excess at energies of 8 to 50 
Bev has also been obtained. The results, while 
in general agreement with those of other work- 
ers, confirm the expected tendency for the pos- 
itive excess to decrease with energy. 


*Now at Stevens Institute of Technology, Hoboken, 
New Jersey. 

tNow at Osaka City University, Osaka, Japan. 

tNow at University of Michigan, Ann Arbor, Michi- 
gan. 


NEUTRAL LEPTON CURRENTS AND NEUTRINO 
DETECTION. R. W. King, Physics Department, 

Purdue University, Lafayette, Indiana (Received 
October 6, 1960). 


The interactions of neutrinos with complex 
nuclei are explored assuming the existence of a 
neutral lepton current. 


ELASTIC SCATTERING OF 146-Mev POLAR- 
IZED PROTONS BY DEUTERONS. Herman 
Postma* and Richard Wilson, Cyclotron Labo- 
ratory, Harvard University, Cambridge, 
Massachusetts (Received October 7, 1960). 


The polarization and differential cross section 
of 146-Mev protons elastically scattered by deu- 
terons have been measured in the center-of-mass 
angular range from 3.9° to 170°. A description 
of the liquid deuterium target and detection ap- 
paratus which permitted the measurement of the 
elastic events over these angles is given. At 
small angles the proton was energy-analyzed; 


at large angles the recoil partner of the proton, 
the deuteron, was identified and energy-ana- 
lyzed. Although the energy resolution of + 1.7 
Mev and the angular resolution of + 2.0° were suf- 
ficient to separate quasi-elastic events from elas- 
tic events at most angles, they were insufficient 
to resolve unambiguously the 2.3-Mev inelastic 
events resulting from the formation of virtual 
deuterons. The measured cross section is in 
qualitative agreement with cross sections at 
neighboring energies; no comparison of meas- 
ured polarizations is possible due to the lack of 
other experiments. The cross sections and 
polarization for angles less than 80° c.m. are 
well fitted by the Kerman, McManus, and 
Thaler theory using the Gammel and Thaler 
nucleon-nucleon potential description. The en- 
ergy and angular dependence of the large-angle 
pickup cross section proposed by Chew and 
Goldberger agree well with experiment; how- 
ever, the small measured positive polarization 
is not predicted by this theory, and probably in- 
dicates destructive interference between the 
direct and the pickup scattering. 


*Now at Oak Ridge National Laboratory, Oak Ridge, 
Tennessee. 


SYMMETRY OF THE ND~' SOLUTIONS FOR 
COUPLED SCATTERING AMPLITUDES. J. D. 
Bjorken, Institute of Theoretical Physics, Stan- 
ford University, Stanford, California and M. 
Nauenberg, Institute for Advanced Study, Prince- 
ton, New Jersey (Received October 11, 1960). 


A proof is given that the ND~* matrix solutions 
for coupled scattering amplitudes are symmetric 
if the discontinuity of the scattering matrix 
across the unphysical cut is symmetric. 
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